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DETERMINATION OF THE MOBILITY PROFILE 
IN OaAs-MESFBTs 

W. Prcat 


1 Introduction 


/!♦ 


1 , 1 Motivation and Purpose 

The Institute for Semiconductor Technology (HLT) at the 
University of Duisburg (West Germany) produces and investigates 
galliura-arsenide-metal-semiconductor-f ield effect-transistors 
(GaAs-MESFET) . It is the goal of the Institute to optimize 
certain characteristics of the structural elements through 
directed influence on technology. 


Technology Physical Parameters Characteristics of 

Structural Elements 



Fig. 1.1.1: Optimum cycle of structural element production 


The producer of the structural elements observes the influence of 
technology on the physical parameters and on the finished /_2 

product. A change in structural element production aimed at 
optimization of the desired relationship requires full awareness 
of all available data. The interdependence of these quantities 
is so great and varied that conscious optimization is 
unthinkable without awareness of the physical parameters. 

The object of this thesis is to supply producers of structural 
elements with a measuring process which incorporates one of these 
physical parameters, namely charge carrier mobility, and its 
local dispersion in the channel (mobility profile). 


♦Numbers in the margin indicate pagination in the original text. 
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It is already possible to determine charge carrier concentration, 
DC power characteristic curves, and microwave characteristics by 
means of measuring technology. The incorporation of the mobility 
profile, therefore, supplements the Institute's measuring 
technology. 

Due to the reciprocal action of all data involved in the 
production process, it is advantageous to be able to measure the 
physical parameters on the structural element itself rather than 
on a special measuring structure (for example, Hal 1-Kreuz ) . 
Further, the measuring process shored be nondestructive, so that 
the structural element is available for further measurements and 
can be characterized in detail. Both of these criteria led to 
the selection of the measuring procedure. 
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No 

. Methods 

Locally 

resolved? 

Destruc- 

tion? 

Special 

Measure 

Structure? 

Literature 

1 

Chemical 

indirect 

yes 

yes 

[7] 

2 

Long-channel 

indirect 

no 

yes 

[5] 

3 

van der Pauw 
(diff. ) 

indirect 

no 

yes 

[6] 

4 

Magneto- 

resistance 

indirect 

no 

no 

[4] 

5 

Magnetotrans- 

conductance 

direct 

no 

no 

[1] 


Table 1.1.2: Survey of procedure aimed at a locally resolved 

mobility profile 


We see from Table (1.1.2) that procedures 4 and 5 satisfy the 
aforementioned assumptions. The decided advantage of No. 5 lies 
in the fact that mobility b(xg) is measured in a very thin layer 
(4=0, 6-6 nm) at xq, while all of the other procedures determine 
b(xp) as integral mobility over the still open channel of xq to 
x=dp, the boundary of the semi-insulating material (SI). 

For this reason, only procedure No. 5, MAGNETO-TRANSCONDUCTANCE- 
MOBILITY-PROFILING, is capable of definitely recognizing the 
extremely important boundary layer effects at the gate contact, 
as well as those on the boundary surface to the SI. 
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MAQNBTO-TkANSCONPUCTANCB-MOBILITY-PROFILINQ 


1.2 


IA 


This concept identifies the measuring process which evaluates the 
change in transconductance in a magnetic field in order to 
determine charge carrier mobility in the channel. This procedure 
does not determine charge carrier mobility as a mean value over 
the entire density of the channel, but rather as a value on a 
defined position xq. 



x: Channel depth coordinates 

x=0 : Gate contact 

x = dg: Boundary to SI 


Fig. 1.2.1; Schematic mobility profile 

The procedure is applicable to all normal MESFETs with a gate 
length in the range of i >im as used in microwave circuits. 

2 Physical basis /5 

2 . 1 Mobility, transconductance 

If an external voltage U is applied to a semiconducting crystal, 
an electric field strength E occurs in the crystal and applies a 
force to the electrically charged particles. 

"F = q ■ ~E (2.1.1) 

If the charge carriers are not connected to the grid, this load 
leads to a movement of the charge carriers and also to a flow of 
current. Charge carriers do not move uniformly. They are 
disturbed in their movement by collisions with grid atoms and 
ionized atoms. By taking the mean, one arrives at the following 
relationship: 
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( 2 . 1 . 2 ) 


v = at R • g 


v .• Drift speed 
t R : Relaxation time 

m* ! Effective mass of the charge carrier 

Drift speed is proportional to the electric field. The 
proportionality factor is called charge carrier mobility, 

ii-R = b (2.1.3) 

m* 


From equations (2.1.2) and (2.1.3) it follows that: 
v = b • ¥ 


(2.1.4) 


Ohmic law is applied to attain the relationship for the /6. 

conductance of the crystal. 

g = a • E l f 

= q • n • v (Ohmircr law) 

a = q * n • b (2.1.5) 

g : Electric current density 

n : Charge carrier concentration 

a : Electric conductance 


The relationship (2.1.2) is not valid for all high field 
strengths. Speed saturation occurs beginning at a defined field 
strength E 4 o . 

v = v Sat (2.1.6) 

For all field strengths E>Eg relationships (2.1.3) to (2.1.5) are 
no longer valid. This is of great importance in determining the 
measuring process, since it follows that charge carrier mobility 
can only be determined in the valid range of the ohmic law. 


-^gs 


Uds 



Fig. 2.1.1: Output characteristic curve field of a MESFET 
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The measurements must, for this reaaoni always be carried [1_ 

out at working points (U{jg, u fls) * which guarantee E<<Eg. 

The more uniform the crystal structure, the greater the mobility 
of the charge carrier in the crystal. Dislocations, 
discontinuities, and other grid disturbances impair the mobility 
of the charge carriers. Also the necessary doping of 
semiconductors is a disturbs the crystal grid and leads to a 
reduction of charge carrier mobility. For this reason a 
compromise has to be found between high current density resulting 
from high discontinuities and high mobility resulting from low 
discontinuity concentration. 

The measuring procedure depends upon the transconductance of the 
magnetic field. For this reason, we shall now concentrate on the 
concept of transconductance. 

S — Livn 


ips 


In the linear range of the characteristic curve field (saturation 
ranges no longer being considered) transconductance increases 
linearly with Uq s . 



Fig. 2.1.2: Output characteristic curve field in the 

starting range 


For each respective working point the transconductance /8_ 

characterizes the grid penetration factor contained in a change 
of the input quantity Ugg upon the output quantity log. The 
transconductance of a MESFET is determined by gate capacity and 
current depth. In the linear range a relationship between 
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transconductance and electric conductivity can be obtained with 
the help of equation (2.1.7) and the ohmic law. 


S . Of • tim 
\ &V 0 ^ °GS 


( 2 . 1 . 8 ) 


E : Field strength in the channel 

A i Cross section of the channel opened through Uq§ 


Transconductance is proportional to a in the linear range. 


2 . 2 Physical basis of the measuring principle 

Aa already mentioned in Chapter (1.2), the measuring principle 
essentially consists of measurement of the tansconductance 


— with magnetic field 

— without magnetic field 


In order to better understand the affects and the derivation of 
the algorithm! which, on the basis of these measurements, 
determines charge carrier mobility, we shall now investigate the 
influence of the magnetic field upon transconductance. Since 
equation (2.1.8) makes a proportionality between transconductance 
and electric conductivity evident, it suffices to describe the 
magnetic field dependency of o. 

In the electric and magnetic field a force F acts upon charge 
carriers in the crystal, which in practice is composed of an n- 
doped material or electrons: 

~F = -e * (E + ~v * B) (2.2.1) /9. 

e : Electric elementary load 

B : Magnetic flow density 

Under the influence of this force, the charge carriers drift away 
from their course, which would run parallel to the electric field 
if it weren't for the magnetic field. On the one side of the 
conducting crystal a relative charge carrier deficiency occurs, 
while on the other side a relative charge carrier excess comes 
into being. A potential difference occurs between the two sides, 
the "Hall-potential" Ujj. 



Fig. 2.2.1: Geometry of a Hall structure 


The field streng_th_J2y increases until it has compensated the 
Lorentz portion v* B . Statically, no force acts upon the charge 
carriers * ■> direction 'ey. They cross the semiconductor along the 
field '.'.Li've of E x , that is the shortest course. Along the 
Current depth has changec^ in opposition to the situation B=0. It 
increases in direction -cTy. No change in electric conductivity 
is to be expected. A reduction of conductance first occurs when 
the transient phenomenon is disturbed between Lorentz /10 

portion VxB and E y and the charge carrier no longer crosses the 
crystal by the shortest course. There are two causes for this: 

— the physical resistance effect P (n,p,b n ,bp) 

— the geometric resistance effect G R (L/W, 0 |j) 

* ( 2 , $ ) ( 2 . 2 . 2 a) 

E : Resulting field strength in the crystal 

p : Concentration of positive charge carriers 

With these assumptions a statement of electric conductivity is 
defined as a function of magnetic field dependency : 

o(B) = o ( 0 ) (2.2.2b) 

G r (L/W, © H ) • P(n,p,p,b n ,b p ) 




n 


1 1 

ll 
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2.2.1 Phyaical conductance reduction 

In the case of a two-unit conductor, n and p are equal or of 
equal magnitude and a compensation of the Hall field strength 
occurs. Positive as well cs negative charge carr^ers^are 
displaced in the same direction by the force q • v * B, so that 
on the sides charge carriers are simultaneously charged 
positively and/or negatively. Compensation occurs and the Hall 
field strength decreases. Since, however, the Lorentz force 
continues without decrease in strength to act upon the charge 
carriers the current lines become crooked. The path through the 
crystal becomes longer and electric conductivity decreases. / 1 1 

Aside from the concentrations n and p, the effect is still 
dependent upon the charge carrier mobility (p. [2]). For n>>p, 
not enough p-charge carriers are available to cause a noticeable 
compensation. For doped semiconduct er material it is, therefore, 
valid, that: 

P(n> >p) = 1 (2.2.3) 

Accordingly, the physical resistance effect cannot be responsible 
for the decrease in transconductance of highly doped MESFETs in a 
magnetic field, 

2.2.2 Geometric conductance reduction 

Every real Hall structure has finite dimensions, which are 
defined here by 

L : Length in direction of control current 
W : Length in direction of the Hall potential which is 
formed 

(see Fig. 2.2.1). Along W the Hall structure is metallically 
contacted, to facilitate better control current feed. Near this 
metallic contact: 

Ey = 0 

On ideally conductive materials, field curves of the electric 
field strength will always be perpendicular. This means, 
however, that no Hall field strength can develop there and that 
the current lines are therefore crooked, /12 
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Fig. 2.2.2: Current and equipotential field curves of a 

semiconductor plate in a magnetic field with 
ratio L/W = 2.4 according to /2/. 

Figure (2.2.2) clarifies how the Hall field strenght disappears 
and causes the current curves to become crooked at the contacts. 
In the middle of the Hall structure the Hall field strength is 
completely developed and the current curves run parallel to the 
applied electric field E x . 

In the case L/W-*“>, the ranges with crooked current curves are 
relatively minimal and the effect can be ignored. 

Gr ( L/W-* ~) = 1 (2.2.4a) 

Actually in an actual short channel-MESFET the relationship is 
approximately: 

L/W = 0.0023 

For this reason we must consider the boundary case L/W = 0 more 
closely. In this case the transient phenomenon is completely 
destroyed and therefore 

/13 

G R ( L/W 0) = Gr max (2.2.4b) 

The mathematics for the description of this process are quite 
involved. Here it will have to suffice to quote statements and 
results ((2] p. 74ff). 


Generally, Maxwell equations with a problem set up according to 
Figure (2.2.2) are applied for universal L/W. Kuhrt and Lippmann 
find a solution through conforming illustrations: 



W _ 1 Z.-3 

L cos © H N 


(2.2.5a) 


Zr and N are integral expressions in accordance with: 



(2.2.5b) 
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(2.2.5c) 



dir 


•/, , - k (VT^) 

N ~ 1 - k'(k) 


( 2 . 2 . 5d) 


K(p) ia the "complete elliptic integral". Equation (2.2.5) 
describes the geometric function G R for universal L/W. For the 
application! however, we must consider the boundary case G H (I/W = 
0 ). 

The boundary value formation l/W-*0 on the equation (2.2.5) 
applies the following geometric function in a simplified form: 

Or (b/W 0 »@h) = 1 + tan 2 @R (2.2.6) 

As the form of the equation (2.2.5) already suggests, the 
boundary value formation is not elementary, but rather requires 
complex calculations. Here, too, one should refer to literature 
/ 2 /. 

The equation (2.2.6), however, is a fairly manageable 
mathematical intermediate result. The variable © u is 
nonetheless a disturbing factor. It still has to be replaced. 

For this purpose one uses the relation: 

tan H = o (B) * H h (B) ■ B (2.2.7) 

This relationship is known from the theory of the Hall 
generators. For strongly doped semiconductor material in the 
exhausted range Rr(B) becomes: 

R H = - _A (2.2.8) 

e * N d 

Np : Doping level 

Kuhrt and Lippmann report, that this is valid for magnetic flow 
densities up to 18 tesla ([2], p. 71). In equation (2.2.7) the 
physical conductance reduction is indicated with a - a(B), which 
has already been excluded for highly doped materials in Chapter 
(2.2.1). With the universal definition equation of the 
conductivity for n-conducting material 

o = -e * b • Np • (2.2.9) 

and equation (2.2.8) equation (2.2.7) becomes: 

tari0H = b • B (2.2.10) 
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The mobility of the charge carriers can be calculated using /15 
equations (2.2.8) and (2.2.9). 


b = R I{ • o (2.2.11) 

It is particularly important here that the Hall constant R R and 
electric conductivity are not functions of magnetic flow density. 
It follows that: 


b # f ( B ) (2.2.12) 

Gquation (2.2.10) inserted into equation (2.2.6) results in: 
G r ( 1/W->0) = 1 + b 2 • B 2 


o(B) 


O 

1 + b 



B* 


(2.2.14) 


With respect to b, equation (2.2.14) can be solved thualy: 



Since b is not a function of magnetic flow density, equation 
(2.2.15) can be considered a physically viable possibility for 
determining charge carrier mobility. Prerequisite to this 
application is the preservation of all idealizations which were 
made in attaining equation (2.2.15): 

1 . n > > p 

2. L/W -♦ 0 

3. The field curves of the electric field strength are 
perpendicular to the contacts. 

Ey = 0 

4. To solve the integral expression (2.2.5) a Taylor series /16 
development of the course of an electron was used. Between 
two collisions with the grid the electron’s course in the 
magnetic field describes a cycloid [2J. The following 
simplifications are therefore valid, if only a small part of 
the cycloid is traversed, so that the approximation using 
Taylor series developments is sufficiently precise. In the 
literature the following boundaries are cited: 

b • B T 1 [4] 

b * B << 1 [2] 
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2.3 


Application of physical phenomena to real 
structural elements 


LXL 


2.3.1 QaAs-short channe 1-MESFET 


Let us now consider a schematic drawing of the function and the 
geometry of the structural element. 



Fig. 2.3.1: Schematic cross section of a MESFET 

(R) = space charge zone 
(SI)= semi-insulation material 


Typical geometric data are: 


Gate length 
Gate width 

Gate-source distance 
Gate-drain distance 
Total channel density 
Passage area height 


L = 0.7 jjm 
W = 300 ,pm 
IqS = 2.5 pm 
Iqq = 2.5 >>m 
dp = 60 to 300 nm 
h = 200 to 400 nm 


Source (S) and drain (D) are block-free ohmic contacts. /18 

Gate (G) is a blocking metal-semiconductor control electrode. 
Current flows from the source to the drain and is directed by the 
gate electrode. To this end voltage U GS is applied between the 
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gate and the source. This voltage determines the density of the 
space charge zone under the gate and with it the density of the 
open channel. It can vary between: 


U T < U GS < U Bi 

U T : the channel is constricted and the voltage applied between 

the drain and source cannot drive any current (except for 
a residual current) through the structural element. 

Ugi : the space charge zone is discharged. The channel is open 

from the gate contact to the SI, so that the channel 
resistance is minimal and the drain-source current I|)g is 
maximal . 


Density and conductance of the open channel range can be 
determined as follows [12]: 

<**«*<.- 

^0 ’ (2.3. ta) 

U ? * U 2 i - Ur (2.3.lb) 


r _ q • n ■ 6 * d ■ W 

b l< " 1 (2.3.2a) 

G = 6 • ~7 — (2.3.2b) 

K L 

d : open channel density 

e q : di-electric constants of the vacuum 

e r : relative di-electric constants from GaAs (11 to 13) 

Up : total potential difference between completely opened /19 

and constricted channel 
N d : doping level in the channel 

2.3.2 Structural element as Hall structure 


According to Figure (2.3.1), the structural element is exposed to 
a homogenous transverse magnetic field. As a result the electric 
conductivity in the entire area between the drain and source is 
influenced. Evaluation of this conductivity variation as 
described in [4] is not limited to the channel, but rather 
includes the contact areas and the areas between the contacts. 

The influence of these dimensions diminishes according to how the 
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passage and contact resistances become negligible with regard to 
channel resistance. At voltages U QS > 0 Volts the channel is 
wide open. Then, in spite of gate lowering, channel resistance 
is dimensionally equal to passage and contact resistance. 

By superimposition of a small alternating voltage Ug a on the grid 
potential of the gate U G c it is, however, possible, not for the 
entire open channel density d, but solely for the channel layer <J 
modulated by Ug B to be considered an effective Hall structure. 

It is necessary only to evaluate the additional alternating 
current ijg flowing through Ug B . 

L 20 


l 


X 


L : Channel length 
W : Channel width 
<5 : Layer density 


Figure. 2.3.2: Tolerance of effective Hall structure 

In preliminary approximations it may be assumed that the length 
and width of the channel concur with those of the gate contact. 
For this reason the values (L = 0.7 W = 300,ym) taken from 

Figure (2.3.1) can be used. The channel layer density is 
calculated with half of equation (2.3.1). 

<1 = d(U GS - u gs ) - d(U QS ) (2.3.4a) 

6 c 


Calculating on the basis of the MESFET data from Figure (2.3.1) 
and using the values: 
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U BT U GS 





N„ = 0.5 • 10 18 cm -3 
u gs = 50 mV 

the density of the channel layer in the proximity of the gate Qg 
= 0.5 V) is found to be: 

6 - 2 nm /2 1 

In the proximity of SI (Uqs = -1.3 V) it is: 

<J = 0.6 nm 


For minimally doped MESFETs the density of the layer can increase 
to 6 nm. In comparison with the density of the right-of-way area 
(80 to 400 nm) these values are minute. Since resistance is 
reversely proportional to density, channel layer resistance is 
much greater than passage area and contact resistances. 

On the basis of the geometry of Figure (2.3.2) a side 
relationship occurs, given as: 


L/W = 0.0023 


This value very nearly equates to the boundary case L/W = 0, 
which was used for the simplified evaluation of the geometric 
conductance reduction. The influence of the magnetic field on 
the conductivity of a thin channel layer in the structural 
element can, therefore, be described sufficiently well with the 
following expression (Bee Chapter 2.4): 


o (B) 


g (°l 

lm 


Y 


(2.3.5) 


2.3.3 Influence of magnetic fields on transconductance 

Up to this point it has been possible to determine the /22 

magnetic field dependency of the electric conductivity. The 
Magneto-transconductance-mobility-profiling procedure does not, 
however, depend upon conductivity measurements, but rather upon 
transconductance measurements. For this reason it is necessary 
to investigate how transconductance changes in magnetic fields. 


According to equation (2.1.7) the definition of transconductance 
is realized by means of superimposition of the alternating 
voltage Ug s on the grid potential of the gate Ugg : 


5 = i Lrv l 


* ^CtS 




(2.3.6) 


Chapter (2.3.2) already introduced the alternating voltage of the 
gate Ug s . It is necessary because it is used to evaluate a thin 
channel layer (locally resolved mobility profile) and the 
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influence of the right-of-way and contact resistances can be 
reduced. A transconductance measurement in accordance with 
equation (2.3.6) has the advantage that it produces a dimension 
which 

— is only dependent upon ij B and Ug a . 

— guarantees that the evaluation limits itself to the thin 
channel layer. 

— is plain and sensible. 

The magnetic field dependency of the transconductance is 
determined by ijg. The working point and the alternating current 
of the gate are independent of the magnetic field. 

ids(t) = ^ K (6, t) U DS (2.3.7) 

In equation (2.3.7) it was taken into consideration that /23 

— the work point lies in the starting range of the MESFET, so 
that the ohmic law is valid. 

— the time independence of the current i^g can be related back to 
the channel layer which has been modulated in the cycle of Ug s . 

If one takes into consideration that the channel conductance is 
dependent upon the electrical conductivity in accordance with the 
following: 

G k (<J) = o • (5ft) * W (2.3.8) 

L 


then the following holds true for the relationship of the 
transconductance without magnetic field S(0) to the 
transconductance with magnetic field S(0): 


SjLOl 

S(B) 


idslOl = £101 
ids<B> °<B) 


(2.3.9) 


The magnetic field dependency of the transconductance concurs 
with the dependency of the electric conductivity determined in 
Chapter (2.2). Integration of equation (2.2.15) leads to the 
measurement equation: 

b = -T (2.3.10) 

B V S(B) 


2 . 4 Influence of idealization 


/24 


Chapter (2.3) showed that the geometry of a GaAs-short channel- 
MESFET is suited for the utilization of the effect of the 
geometric conductance reduction. The effect of the conductance 
reduction could be carried over to the transconductance 
measurements. Cognizant of the geometry and physics of the real 
structural element, we shall now make a qualitative and, if 
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possible, quantitative analysis of errors which occur when 
adapting a MESFET to the idealizations of equation (2.2.15) 
and/or (2.3.10). 

2.4.1 Finite geometry 

Figure (2.3.2) presents a portion of the structural element which 
is responsible for the magnetic field dependency of the current 
id s . We have established the side ratio: 

L/W = 0.0023 (2.4.1) 

The equations 

GR (L/W, b , B) * 1 + b 2 B 2 (2.4.2) 


° ( 0 ) 


o (0) 

1 + b 2 B 


t: 


(2.4.3) 


are valid under the given: 
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L/W~*0 


(2.4.4) 


Equation (2.4.1) presents an error, whose influence on equation 
(2.4.3) must be investigated. The universal solution of the 
geometric function is: 

i *s. 





cos & 


H 


N 


(2.4.5) 


Z R and N are non-analytic soluble integral expressions. For 
reason integration of equation (2.4.1) into equation (2.4.5) 
not lead to an analytical calculation of the mistake 
finite side ratio L/W. 


this 
doeB 
through the 


The following graphic solution of the geometric function serves 
for a qualitative analysis. 
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Fig, 2.4.1; Geometric function for resistance increase of 
right-angled semiconductor plates [2] 

The representation of the geometric function shows very /26 

clearly how near the graph for L/W = 0.0023 lies on the graph of 
L/W = 0. It should, however, be pointed out, that the product 
b • B does not exceed the value 0.5. This value range is only 
insufficiently represented in Figure (2.4,1). The integral 
expressions in equation (2.2.5) are constant, so that the 
statement encountered in the aforementioned figure for a 
relatively high b • B, can be carried over to a lower value of 
this product. For this reason a noticeable error caused by the 
real side ratio L/W on equation (2.4.3) can be avoided. 


2.4.2 


Passage and contact resistance 


For short-channel MESFETs in starting range the ohmic law is 
valid between the applied voltage Ups and its resulting current 
IpS* The drain-source current may be characterized by the 
following representation of a resistance replacement switch. 
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Fig. 2.4.2: Replacement switch representation characterizing 

the drain-source current 

Channel layer resistance modulated by Ug a /27 

Open channel resistance determined by U Q g 
Contact resistance of the drai.n/source contact 
Passage resistance between drain contact and channel or 
source contact and channel 

The channel conductance in Figure (2.4.2) is represented by the 
parallel connection of Rj((tf) with Rj((d). The superimposition of 
the gate’s grid potential Uq S and its alternating voltage Ug s 
must be considered here. Since only the current portion ijg is 
drawn upon for measuring transconductance, channel resistance 
Rj{(d) remains without influence. The following simplified 
representation of a replacement switch characterizing the 
alternating current ij s results if the contact and passage 
resistances of the drain and source are combined into one contact 
and passage resistance. 


S D 




Fig. 2.4.3. Simplified representation of a replacement 
switch characterizing ijg 


On the basis of Figure (2.4.3) the influence of errors caused 
passage and contact resistance can be investigated. The 
following relationship serves as a solution statement: 


SfB? ft Co) 

SCO) "W 

S LI) _ If +■ ftfl f K<*) 

S(£>) " l + A.U+6‘8 1 ; ♦ A^G+tJa 1 ) 

£ 


(2.4.6) 


by 
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bg ! Charge carrier mobility in the channel /28 

bg ‘ Charge carrier mobility in the paeaage area 

b c : Charge carrier mobility in the contact area 

This statement begins with the assumption that every resistance 
is subject to a geometric c6nductance reduction in the boundary 
case L/W =0. A side ratio of 0.01 is not exceeded for the 
passage areas, so that the approximation is sufficiently 
fullfilled in all cases. Taking the contact area as a Hall 
structure is more problematic. Here, current flow and magnetic 
field are not always perpendicular to each other. It is also 
unusual to apply Hall structures to ohmic contacts. 

Nevertheless, the electric conductivity, which is described bb 
analogous to semiconductor materials by means of mobility and 
concentration of the charge carriers, is completely dependent 
upon the magnetic field. The resulting side ratio concurs 
approximately with that of the passage areas. 


For consideration of errors, therefore, a statement based on 
equation (2.6.6) is legitimate and can describe the influence of 
parasitic resistances on the equation 

S(B), 1 . . (2.4.7) 

S(0) 1 + b^B 2 


An evaluation of equation (2.4.7) delivers an uncorrected value b 
for mobility in the channel, while an evaluation of equation 
(2.4.6) delivers a corrected value, the actual mobility in the 
channel. The error which is being sought results from an 
evaluation of equation (2.4.7). To this end the two equations 
are equated to each other: 

_J _ jjc. + <3 

' 1'U H c l S l ) *■ ttyUMW 

Transformations of the equation result in: 


(2.4.8) 
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(2.4.9) 


A simplification for small k can be carried out by means of 
Taylor series development. It is based on: 

VTTTT 1 + 0,5-k 


This approximation car, be 
geometry guarantees, that 


applied to equation 

R k ( 6) » Rg«R c 


(2.4.9) , 


since the 



- it i - Vc - -h- 

1 t S..(6) y <*/ i 


K <6 >' 


2a_ 

w 


•fr- 


fa. 

6 l 


(2.4.10) 
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The error is defined at: 


f 


b 



1 


A comparison with equation (2.4.10) gives us: 




± 

l 



Discussion of equation (2.4.12): 


(2.4.11) 


(2.4.12) 


Equation (2.4.12) consists of two correction terms. The left 
describes the influence of the contact areas, the right those of 
the passage areas. The direction of the influences is given by 
the relationship of the mobilities. The resistance relationship 
only has influence on the strength. 


Ohmic contacts generally have a high concentration and a low 
mobility of charge carriers. bg is therefore generally smaller 
than b, so that the first term is negative. 

In passage areas the mobility is higher because of the /30 

larger distance from the surface, so that the second term will be 
positive. The errors do not generally add up. 

The resistance Rg(d) is a function of the depth in the channel. 

At the gate contact <5 is large and therefore resistance ia small. 
As distance from the gate increases, 4 decreases and resistance 
increases. For this reason the influence of mistakes caused by 
resistance relationships is greatest at the gate contact. 


In the following we shall assume a quantitative mistake 
estimation. If we limit ourselves here to the passage areas, 
then this will be an estimation upwards, since the errors from 
the contact areas must be substrf ted quantitively from it. 


f b< 



(2.4.13) 


At the gate contact b < bg, Depending upon the profile, a drop 
in voltage up to half of the maximum value can result. The 
mobility bg is an average dimension above the height of the 
passage areas and is not dependent upon the working point. From 
the primary approximation the following can be set: 

b max ~ b B 


The maximum mobility in the channel corresponds to the mobility 
of the passage areas. With the aforementioned assumptions, then, 
it is valid that: 
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The 


b min = °- fi b B 

resistances are calculated according to equation (2.3.2): 


1 B 

h 


R k (6) 


1 

n • b • q 


L 

W-6 


1 

n-bg.q 


2>i B 

W-h 


(2.4.14a) 

(2.4.14b) 


Length of the passage areas between contact and /31 

channel 

Height of the passage areas 


The factor 2 in the equation of the passage resistance considers 
the addition of the passage resistance between the gate and 
source and the gate and drain. Integration of the values 
L = 0.7 .ym and a layer depth 6 = 2.6 nu results in a resistance 
relationship: 


3 _ j b ' & ' 

V*) ’ ^ L 


o, mb 


(2.4.15) 


Integration in equation (2.4.13) gives us an upward estimation of 
the mistake at the gate contact at 

f 0 = -14% 


This error decreases through the influence of the contact area. 

As distance from the gate increases, the number of errors 
decreases rapidly, since 

— layer density d decreases. 

-- b increases. 

In the proximity of the semi-insulation material b decreases more 
drastically again, so that errors, in spite of higher resistance, 
Increase again. 

Exact error estimation is only possible on the basis of 
corresponding structural element data and their distribution in 
the channel. The principle influence will be examined in the 
following. 
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Fig. 2.4.4: Distribution of the error function f^ and its 

causes in the channel 

An error boundary is cited in literature source [1]. This value 
is maintained in wide portions of the profile. A value of 10* 
may not be exceeded in the boundary layers. 







2.4.3 Idealization of the geometry and summary obaervation of 

physical and geometric error Influences /33 

In Chapter (2.3.1) the area under the gate was defined as an 
active Hall structure. The geometric conductance reduction for 
this structure was then calculated according to 

SfB) . iiSl 

M'B 1 

With the exception of the approximation L/W~*0, the Hall 
structure must have metallic contacts in order to fulfill the 
aforementioned equation ideally. A short circuit of the Hall 
structure and/or a forced right angle position of the electric 
field strength on the electrodes can only occur by means of very 
high contact conductivity. Establishment of such an ideal 
contact is not given in the MESFET. Nevertheless the connection 
of the thick passage areas (hs^200 nm) onto the thin channel 
layer (<J**2.6 nm) can be considered a relatively ideal contact. 

In the passage areas the electric current density is very much 
smaller and therefore 

E yB = r H * «xB ‘ ° . 

The y-components of the electric field strength are very small. 
Further, the geometric resistance effect weakens E y in the 
passage areas with 

1/W = 0.008 

The connection from the passage area to the channel can be looked 
at as a boundary layer free of space charge, so that it is true 
that : 

divlT = 0 (2.4.16) 

The tangential components of E are constant. Since the y~ /34 
components of E are near zero in the passage area, the y- 
coraponents in the boundary layer are minute as well. 

An even more thorough consideration of these and other physical 
sources of errors would be beyond the scope of the present study. 
On the basis of the observations we have made thus far it is to 
be expected that according to 


4 s - iT'I/m) “ 1 (2.4.17) 

electron mobility can be measured with a high amount of 
precision. Nevertheless this relationship can only be valid with 
a high amount of surety if the experiment confirms its resultant 
equation : 
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S(B) * 


(2.4.18) 


5(0) 


The result of such an experiment is shown in the following 
figure. 



Fig. 2.4.5: Experimental confirmation of the evaluation 

equation 

Over the magnetic flow density B we add - i ' . A comparison 
with equation (2.4.17) results in a line equation with the 
gradient b. 



(2.4.19) 


If the graph in Figure (2.4.5) is an even one, then it is a proof 
for 


— the quadratic relationship between S(0) and S(B). 

— the independence of the mobility b from B in the exhausted 
range of the magnetic field. 

Experimental control of the resultant equation is needed because 
of the varied causes of errors. An exact calculable /36 

control of the error function fj-, according to equation (2.4.12) 
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is not useful because of the high technical effort required to 
attain the structural element data which depend on the test and 
each respective working point. On the basis of the foregoing 
considerations the error boundary caused by physical idealization 
can be established at +/-5X. 


Measuring room 


Measuring principle 


In the following chapter we shall first summarize the resultant 
equations for the measuring procedure and then discuss the set up 
and the functioning of the measuring room. 


To determine mobility in the channel we have the resultant 
equation 

(3.1.1) 



The equation determines mobility at a depth determined by the 
grid potential of the gate Ugg . By superimposition of the 
alternating voltage Ug s on the grid potential Uq S the thickness 
of the space charge zone dpj is modulated to its mean value. 

There is an alternating current ijg which originally is connected 
with the modulated channel layer in the boundary between the 
space charge zone and the open channel. The transconductance 
measurement 



(3.1.2) 


is strongly related to the modulated channel layer. A second (or 
more) transconductance measurement results when the test is 
subject to an homogenous magnetic field of the flow density B. 
With the data 


— S ( 0 } 

— S ( B) 

— B 


equation (3.1.1) yields the mobility: 


b = b(x 0 ) (3.1.3) 

With the agreements of Figure (2.3.1) it holds for the /38 

depth coordinates xq that: 

x 0 = (3.1.4) 


By varying the work points 
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(3.1.5a) 

(3.1.5b) 


0 < X < dg 


u Bi > Uq S > u T 

one can in this manner measure mobility at every depth of the 
channel. For clarity the measuring process is shown in the 


following flow diagram. 



Key: 1) Start 2) Setting of 
U BS 3) Setting of U QS 

4) Measurement of S(0) 

5) Magnetic field "IN” 

6) Measurement of S(B) 

7) Magnetic field "OUT” 
Evaluation: b(x B ) 8) New 
depth? 9) yes 10) no 
11) End 


Figure 3.1.1: Flow diagram of the measuring cycle 
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The following measuring device is required to carry out /40 

this measuring cycle: 



Figure 3.1.2: Model of a measuring device 


The measuring device consists of: 

- Adjustable direct voltage sources for the work points U GS , 
U DS 

- Alternating voltage source Ug 8 

- Voltage meter for determining alternating voltage 

- A device for neutralizing and determining channel current 

i dB 

To reduce errors the alternating current Ug s in accordance with 
equation (2.4.12) is set at approximately 

Ug a = 50 mV (3.1.6) 

In the starting range of a MESFET transconductance has a value of 
about : 

S = 1 mS (3.1.7) 

Using equation (3.1.2) these two values produce a current: 

Z41 

id s = 50 v A (3.1.8) 

Such a small current can no longer be measured accurately enough 
with a conventional ampermeter to permit its use in the measuring 
procedure. Further, a special filtering technique is required in 
order to separate the current i^g from the equal portion Ipg and 
to reduse the noise effects. For this reason, a lock-in 
amplifier is installed in the measuring room. 
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The technical realization of the measuring device has the 
following configuration: 



Figure 3.1.3: Technical realization of the measuring device 


Application of a load resistance becomes necessary, since the 
lock-in amplifier can measure only voltages. Over the load 
resistance Rj, falls away voltage 

= ids ’ (3.1.9) 

The setting of R^ is so small that the current i^g is not 
diminished. /42 


3 . 2 Measuring equipment and voltage supply 

3.2.1 Lock-in amplifier 

The "Dynatrac 393" lock-in amplifier is the central measuring 
instrument of the measuring device. An evaluation of the 
physical measuring principle is possible only using this type of 
measuring technique. The following figure shows a detail of the 
application range of this measuring instrument. 


35 



0 Lock-in Verstarker 


U 


z) 

Meflteil 

Sinus - 

*ORef. ^Signal 

gen. 

1 

in in 

T 1 


s (t) 
R 


DS 4 U GS 




s (t) 
a 


7; 

Rauschen 


'Ausgang 

(DC) 


Figure 3.2.1: Application range of the lock-in amplifier 

Key: 1) Lock-in amplifier 2) Sine generator 
3) Measuring part 4) Reference in 5) Signal 
in 6) Exit 7) Noises 


In the measuring instrument a comparison between the reference 
signal and the measuring signal takes place according to /43 
the following principle: 

The phase regulating cycle in the lock-in shift phases the 
reference signal Sjj(t) until the product 

s R (t-Tr)«s a (t) 

reaches a maximum. At the maximum state of the product both 
signals have the following per unit time course: 



Figure 3.2.2: Comparison of reference and entrance voltage 

of the lock-in amplifier 
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Over a full period the product totals: 


In the process, 


all portions fall 


i 


T 


out 


I '=*</*•'*> 

t 


*0 1/L) el Jjl 


— whose mean value over T is different from zero (equal portion). 
— whose frequency is a multiple of the basic frequency (harmonic 
distortions) . 


Further, the signal/noise power ratio is optimum, since, 
analogous to the technique of the appropriate filter, signal 
portions with the highest correlation to the reference voltage 
demonstrate a maximum upward integral total. Ideally, white 
noise will be eliminated. /44 

The output of the lock-in amplifier is a normalized direct 
voltage of the form: 

O'R (I-?) ‘/Ogll') j 

(t) | (3.2.1) 



The signal dimension s a (t) is composed of the equal portion Ipg 
R^ and the alternating portion ijg Rj,. In practical use there 
occur alternating currents of 

3 to 300 yA 


which must be measured with an accuracy of approximately 1*. 
Further, the channel resistance Rg(d) lies in the range of 
approximately 30 to 100 ohms. The author is unaware of any 
conventional measuring instrument which is capable of this 
measuring task under the given circumstances. 


The lock-in amplifier cannot measure direct currents, which is 
the reason that a load resistance R^ was applied in the drain- 
source cycle. Analogous to equation (2.4.12) R^ brings about an 
additional error at; 


f. 



' L . JL- 

i 'V*; 


(3.2.2) 


R l is kept at such a dimension (1 ohm), that the error does not 
generally go over 2%. The current ijg can be produced from 





(3.2.3) 


with sufficient accuracy . is the decrease in alternating 

voltage, which causes i<j s on the load resistance R^. 


37 


The exit of' the lock-in amplifier is a direct voltage 
proportional to 


745 


U 


L 


I 


40 V 


(3.2.4) 


Sg is the entrance sensitivity set on the lock-in. It covers the 
value ranges: 


0.1 < S E < 3 V 

The measurement accuracy is given by the producer at: 

1* Amplifier errors 

0.05 %/K Temperature drift 


The filter selectivity amounts to 12 dB/octave. 


3,2.2 Keithley digital multimeter 

t 

The measuring room contains two Keithley digital multimeters, the 
model ”177 Microvolt DMM " and ”179 TRMS DMM” . Both models are 
compatible with IEC-Bus. 


Model 177: 


Assignment : 

Error: 
Entrance : 


Measuring alternating voltage of gate Ug 0 and 
IEC-Bus cut-off position 
0.9% (at 50 mV) 

10® ohms/ 75 pF 


Model 179 TRMS: 


746 


Assignment : 

Error: 
Entrance : 


Measuring exit direct voltage of lock-in and 
IEC-Bus cut-off position 
0.05% (at 10 Volt DC) 

10® ohms/ 75 pF 


3.2.3 Voltage supply 


According to Figure (3.1.2), the direct voltages Ugg and Upg are 
needed along with the alternating voltage of the gate Ug s . The 
alternating voltage generates a sine oscillator in the Dynatrac 
393 lock-in amplifier. 


Sine oscillator: 


Frequency: 
Exit : 

Error : 


0.1 to 10 KHz (here 1 KHz) 

0 to 1 V at 1000 ohms 

2% harmonic distortion, drift (see Chapter 4.1.4) 
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The IEC-Bub compatible two-channel voltage source "Rhode & 

Schwarz NGPS" serves as the direct voltage supply. The following 
data are valid for the "LOW" range. 

Direct voltage source: 

Voltage: +/- 16.835 V 

Current: 100 mA (130 oA current limit) 

Error: So small that the direct voltage supply can be 

considered ideal 

3.2.4 Test circuit /47 

The technical measuring circuit calls for a load resistance Rj, in 
the drain-source circuit and a switch for the feed-in of the gate 
voltage Uq S + Ug s . Since load resistance R^ (1 ohm) is very 
small, already small overvoltages lead to high currents and with 
them to destruction of the structural element. The cause of such 
voltages are: 

manual errors 

- switching the voltage sources on and off 

The following construction was developed for the protection of 
the test and as technical measuring wiring: 



The R-C combination Ry Cg provides for feed of gate voltage Uqs + 
Ug s . Ry prevents a short circuiting of the alternating voltage 
Ug s over the direct voltage source Ug S . The condensor separates 
the direct voltage from the alternating voltage generator. An R- 
Zener breakdown is installed at the gate entrance Uq S . /48 
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Figure 3.2.4: Acceptable value range of the gate-source 

voltage 


The circuit causes the voltage Ua established from U 7 to to be 
equal to the actual voltage at the gate. Voltages which lie 
above this value range are short circuited. The measuring range 
of Uqs looks like this: 

+0.5 V < U GS < -5 V (3.2.5) 

This condition is satisfied when Up = 0.6 V and Ujr = 5.6 V. 

A further protection precaution is the relay Kj. It is turned on 
directly by the measuring program and first establishes the /49 
voltages at the drain and gate when the start-up process has been 
completed. After every measurement the relay once again 
separates the test from the voltage supply. 

These precautions are necessary and effective to protect the test 
against voltage peaks, which occur during switching, particularly 
when turning off the direct voltage source, and can lead to 
destruction of the test. 
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The test circuit is the distributor position for the entire 
measuring construction. It is essentially responsible for the 
ease in inspecting the measuring room. This is mode clear by the 
following figure; 
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Figure 3.2.5; The test circuit in the block circuit of the 
measuring device 

Key according to number: 1) Test wiring 2) 

Voltage source 3) Entrance 4) Exit 5) Lock- 
in amplifier 6) Test 7) Control voltage 8) 
Relay box 

3 . 3 Magnetic field regulating circuit /50 

Figure (3.1.3) shows the most detailed description of the 
technical realization of the measuring device to date. 
Nevertheless, the central calculation control (Chapter 3.4) and 
the adjustment of the magnetic field are lacking. 

It is, however, precisely the adjustment of the magnetic field 
which is of particular importance in constructing the measuring 
room. The exact measurement of magnetic flow density is quite 
involved. Furthermore, there are certain requirements which must 
be maintained: 

— Adjustment of a variety of magnetic field sizes during a 
measuring cycle 

— The most exact adjustment possible 
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— Fully automatic performance of the adjustment procedure by 
means of a terminal 

In Chapter (2.4) we called for technical measuring control of the 
validity cf the initial equation 

5 ( 0 ) 

A + l>'B l (3.3.1) 



This control should take place through verification of the linear 
connections 



(3.3.2) 


m-X = y 


(3.3.3) 


In order to obtain the settings of both magnetic field //jl. 

values (B = 0, B=Bq), which suffice for evaluation of equation 
(3,3.1), many (in practice n=5 to 10) varied magnetic field 
values are set and used to measure the corresponding 
transconductances StB^). 

After the completion of a cycle, n-value groups of the form: 

k ! B k j S ( Bjj ) 

are available. There data are evaluated by means of a linear 
regression (see Chapter 4. 1.3. 2) to a line equation according to 
equat ion (3,3.2). 

This cycle must be repeated for every point on the mobility 
profile. For a complete profile the result is between 50 and 200 
magnetic field adjustments, which by manual adjustment would 
require an operator to be present in the measuring room 
throughout the entire measuring time.* Furthermore, such 
adjustment is very tiring. The error boundary acceptable for the 
adjustment value is subjective and different for every 
adjus tment . 

For this reason the following "digital” magnetic field regulation 
is set for the measuring device. 
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Figure 3.3.1: Block circuit for digital magnetic field 

regulation 

Key: 1) Output voltage source 2) Bl Hall 

probe 3) Bus coupler 4) CBM terminal 

3.3.1 Magnetic- field regulating circuit instruments 

The regulating circuit consists of three groups: 


Group A; Generating elements 

Al: Bruker magnet B-E 10B8 

A2 : Systron Donner Power Supply M7060-30 

Group B: Measuring instruments 

Bl: Hall probe 

B2: Tesla meter Systron Donner 3102 A 

Group C: Adjustment parts and regulators 

Cl: Commander Systron Donner PI 

C 2 : Instrument coupler ICS 4880 

C3 : Terminal CBM 4032 


/5 3 


The elements of group A produce a relatively homogenous magnetic 
field up to a maximum of Bjj, aJ{ = 1 tesla in an air gap of 
approximately 6 cm width. Group B facilitates measuring the 
actual value of B. It consists of a Hall probe, which transmits 
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a magnetic field proportional Hall probe to the second instrument 
of the group, the tesla meter. This instrument iB calibrated on 
the Hall structure and converts the Hall voltage into a value of 
the magnetic flow depth. Magnetic field measurements are 
generally much more difficult and imprecise than electric 
measurements. Nonetheless, the combination used in this 
measuring room are considered quite good. The practical 
precision with which the line equation (3.3.2) can be satisfied 
has helped reduce the justifible prejudices against magnetic 
field measurements. For this reason, the error boundary of the 
producer can be incorporated: 


f T = 1* (3.3.4) 

Group C consists of converters Cl, C2 and the logical element of 
a computer terminal. Cl converts a decimal number from 0 to 9999 
into a control voltage Ug^. of 0 to 10V. U st controls the voltage 
source A2 from 0 to 100% of its maximum voltage. The smallest 
magnetic field variation amounts to 0.01% corresponding to 0.1 
mT. This ensures that the adjustment of every B-value is 
"ideally" possible. C2 allows the IEC-Bus terminal to read the 
digital tesla meter results. /54 


3.3.2 Magnetic field regulation 

The elements of Group A, B and C are connected 
circuit. 


to a regulating 

2 : Storprossen 



Figure 3.3.2: Magnetic field regulating circuit 

Key: 1) Calculator 2) Disturbance dimensions 


The figure shows a technical regulating representation as the 
schematic form of a single regulating circuit. The assignment 
and function of the terminal, the analogous/digital converter 
(A/D converter) and the D/A converter especially are strongly 
schematicized. Blocks A, B, Cg and C R correspond to the 
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component groups listed in Chapter (3.3.1) with the conversion 
factors : 

f A' f B’ ^C c ‘ f C D 


The dimensions with (#) are digital. Considering the great 
effort required for a technical regulating consideration of the 
"dead periods" using data recording and A/D and/or D/A /55 

conversion, dynamic regulation is abandoned. Since A/D and D/A 
conversion can be considered ideal, the problem of conversion 
will not be dealt with henceforth. 

There are two steps in the construction of a regulator circuit in 
accordance with Figure (3.3.2) 

1. Assumption of transfer factors f A and fg, which are given in 
the process and are immutable 

2. Calculation and realization of the transfer factors ^Cc an d 
f r according to the points given in 1. 

C R 

Regarding 1 . ) : 

The transfer function fg describes the technical measuring 
composition and conversion of the actual value of B. Since no 
problems with the A/D conversion or technical measuring errors 
are observed here: 

f B = 1 (3.3.6) 

Transfer function f A describes the conversion of the dimension 
Ugf. into an actual value of B : 

B = f A (U st ) 

Function f A was assumed by technical measurement in /56 

the measuring room: 
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Control voltage 

Figure 3.3.3: Transfer function of instrument group A 

Regarding 2 . ) : 

The processes on the regulator circuit are differentiated 
according to times: 

t < tj^ : Control 
t > t^ : Regulation 

t x : Conversion time of a real value of B in the first, error- 

ridden actual value of B 

In accordance with Figure (3.2.2) the following relationship is 
valid for t < t^: 

Z57 

8 = B S' f C ‘ f A + 1 (3.3.7) 

The term Z is an additional disturbance dimension, which must be 
stabilized. The pure control already works optimally, whenever 
the following is valid: 


B — Bg + Z 


(3.3.8) 


r 
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By means of 
attained: 


comparison the dimensioning specification 



< f A> 


-1 


is 

(3.3.9) 


Transfer function f^ as well as its inversion are not analytical 
functions. For this reason a solution of equation (3.3.9) was 
found by integrating pairs of value from Figure (3.3.3) into the 
statement : 


f c (B) = a^* B + a 2 ‘ B 


(3.3.10) 


The results are: 


a^ = 8.606 • 

10“ 3 V 

(mT) ^ 

a2 = 1.394 * 

10" 6 V 

(mT)“ 2 


(3.3.11a) 

(3.3.11b) 


At t ~ tj_ the real value of the magnetic flow depth was converted 
into an actual value: 


B = B s + Z (3.3.12) 

The switch in Figure (3.2.2) is reversed; the regulator circuit 
is closed. The assignment of part C R takes the dimension placed 
at the entrance: 


AB = B s ~ B 

and generates a correction dimension from it: 

V “* iU st 


/5 8 


The following condition must be satisfied: 

f A (U St + iU St> ' B + 158 

Function f^ is linearized at the working point: 

Q ^ 

~ 0 (3.3.13) 

Now the observations from equation (3.3.7) to (3.3.11) are 
applied to the changes 


* 



B, U st : 


-2k 

0T> 


-l 


(3.3.14) 
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Integration into equation (3.3.10) results in 



&8 • (a | + 2 & 2 ' B ) 


(3.3.15) 


with a^ and 33 in accordance with equation (3.3.11). 

With that the regulator circuit is completely dimensioned. 
Control and regulator functions f c are realized by placing the 
functions into the Measuring program of the computer in 
accordance with equations (2.2. 1C) and (3.3.10). 


Regulation is considered complete when regulator error 


A B “ Bg “ B 

falls below a certain level (1 to 3*). The following figure 
shows the per unit duration of the dimensions affected by 
regulation: 



Figure 3.3.4: Time plan of the magnetic field regulation 


Wi : Waiting time for adjustment of the magnetic field value 

W 2 ! Assumption of the measured value of B 
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The other waiting times (A/D conversion) can be omitted because 
of Wi and W 2 • 


3 . 4 Measuring room automatio n /60 

In the previous chapter a portion of measuring room automation, 
the magnetic field regulation, was described in detail. As we 
saw there, the following points indicate the need for automation. 

— The amount of operator time is considerably reduced, since 
plotting of a highly solvable mobility profile can take many 
hours . 


— Work in the measuring room is very monotonous. This leads to 
a reduction in attention and to error conditions. 


Data recording, data assurance, and evaluation are assumed by 
a computer. 

The ability to reproduce measurement results is considerably 
increased . 


The two last points increase measuring room measurement precision 
and dependability. The computer can calculate even monotonous 
mathematical algorithms, leading to a reduction in errors which 
is absolutely necessary for the measuring procedure (Chapter 
4.1.3). 


The measuring program con be structured as follows: 


m. 



Without operator 


Figure 3.4.1: Structure of computer controlled measuring 

program 


Upon data input, all variable parameters on the special test are 
determined. After completion of this phase, the terminal should 
have received all necessary information and the measuring room 
should be fully prepared to allow the second phase to proceed 
without an operator. An interruption of the measuring program 
should occur only when 


— Boundary data is surpassed 

Errors occur in the measuring construction 

— The adjustment of measuring instruments must be corrected 
manually 
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The last point refers to the limitations of automation. In the 
case of the measuring device, this situation occurs upon 
correction of the entrance sensitivity of the lock-in amplifier. 
If transconductance in the channel radically decreases to SI, 
then the entrance sensitivity of the measuring instrument must be 
reset. Complete automation is possible only if all instruments 
utilized in the measuring room can be controlled centrally from 
the terminal. A programmable lock-in amplifier does not exist at 
the present time. 

Development of automation in the measuring room for magneto- 
transconductance-mobility-profiling has progressed in / 62 

proportion to the availability of measuring instruments which may 
be adapted to the IEC-Bus. In developing programs, particular 
emphasis has been placed on avoiding a complete interruption of 
measurement, with the terminal containing directions for 
resumption of the measuring procedure. 

Example: Initial voltage of the lock-in falls below 20 % of 

total oscillation! 

The technical measuring error threatens to become too 
large . 

1. The computer stops the measurement just before taking in a 
new mobility number in order to avoid errors caused by 
measuring range shifts during an evaluation. 

2. Computer instruction: "Correction of lock-in entrance- 

sensitivity ! " 

3. The operator corrects the measuring range. 

4. Computer asks: "Sensitivity >jV?" 

Operator inputs corrected value. 

5. Measuring program continues! 

This small example may clarify the procedure of logical measuring 
room supervision . We now turn to the complete data technical 
cabling of the measuring room. 
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Key to Figure 3.4.2: 1) Hall probe 2) Test 

3) Magnet 4) Output voltage source 5) 
Voltage source 6) Test circuit 7) Relay 
actuator 8) Lock-in amplifier 9) In 
Reference 10) Exit 11) CBM terminal 12) CBM 
CBM diskette drive assembly 


4 Measuring technique /64 

4 . 1 Brror observation 


In our introduction to the measuring room in Chapter (3) it 
already became evident from the choice o.f modern and highly 
precise measuring instruments that the measuring procedure is 
highly sensitive to errors. Here we deal primarily with the 
imprecisions of measuring instruments and their influence on 
errors in the results investigated. An essential difficulty in 
constructing the measuring room consists of minimizing errors 
caused by the measuring instruments. This task has been solved 
on an engineering level, that is, the dominant error influences 
here have been reduced to a level equal to those of other errors. 

4.1.1 Sources of errors 


Basically 
may cause 
of errors 
and their 
equations 


all of the instruments involved in the measuring room 
errors. In order to attain more meaningful observation 
we direct our attention to the essential parameters 
errors. These paraueters are contained in the initial 


S(B) 


5 ( 0 ) 
4 + 


(4.1.1a) 


5 = Ml 


^ ' $<= Z65 

t L • (4.1.1b) 


Both of these equations result in the evaluation equation of the 
measuring procedure which occurs independently of the error-laden 
dimensions occurring in the measuring room: 



(4.1.2) 


The goal of error observation is to minimize the error boundary 
of b, and, if possible, to record it in numerical values. 
Working from equation (4.1.2) we shall limit ourselves in the 
following to error observation in: 


52 



U^: Lock-in initial voltage 

B: Magnetic flow depth 

Sources of error for these parameters are: 


— Lock-in amplifier 

Hall probe/Tesla meter 


The error influence of the voltage sources for Uqc and Upg can be 
ignored. The same is true for all digital operations of the 
computer and the converters. Because of reference voltage drift 
these will be considered by themselves in Chapter (4.1.4). 

4.1.2 Error sensitivity /66 


Up to the first term, a Taylor succession serves as a starting 
point for the influence of parameter errors on end results: 


B- 


03 


^ Ll0 ^du L (o) f 


o±_ 

ow) 


(4.1.3) 


The partial derivatives can be obtained with equation (4.1.2). 
Here, for relative error: 


b-b 0 

K 


* - 1 
• 4 Ujo) . 

£ UJo)-U l ib) 

1 Ul<°) , 


3 

al/ t (6) 

0 L to) 

U L <*) 


This equation has the form: 




f 


i 


(4.1.4) 


(4.1.5) 


E : Sensitivity factor 

f^ : Relative parameter error 

Relative parameter error f^ corresponds to measuring instrument 
tolerance, while gives the influence of parameter error on the 
error in the end result. Since all of the sum terms can have 
different signs (ideally being normally distributed), the amount 
of the error can be assumed. 
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Vj/lXf,) 2 (4.1.6) 

Under these circumstances an error estimation of f b results, in 
which it is shown that error sources cause a total error of 

f b ^ 0.15 (4.1.6) 

Equation (4.1.7) can be satisfied with the smallest effort, when: 

= | E 2 - f 2 1 . |Ej.f 3 | 

E.'f, = - = 0,0866 (4.1.8) 

All error factors in equation (4.1.4) must satisfy this 
condition. For the term which describes the influence of the 
magnetic field, this boundary, affected by sensitivity factor 
Ei=l, is bound to the condition that the error of the tesla meter 
and the Hall probe not exceed 8.66X, Th<i measuring instrument 
producer lists the tolerance as IX. 

The second term of the equation (4.1.4) is described by: 

r _ 

' l U L (o) - U l (8) 

r __ zUlJo) 

I l w 


The error f ‘2 is dependent upon the amplitude of the lock-in. The 
measuring program guarantees that the lock-in be driven in a 
suitable measuring range and that 


f 2 = 0.01 


(4.1.9) 


not be exceeded. Sensitivity factor E 2 depends on magnetic /68 
field. If f 2 =0.01, an unbalanced equation can be set for E 2 : 



(4.1.10) 


The unbalanced equation (4.1.10) can be satisfied only at 
sufficiently high values of the magnetic flow depth. Inversely, 
an adjustment specification can be derived with this unbalanced 
equation . 
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The following model calculation serves to this end: 


b = 4000 cra^/Vs 
M 0 >/ u l( b ) = i ' 1 

Solution of this evaluation equation according to B and 
integration of the unbalanced equation (4.1.10) results in: 


(4.1.11a) 


(4.1.11b) 


A (VB) 

U l (B) 

f2 and also leads to the 
unbalanced equation (4.1.10) because of the same sensitivity 
factor. The different signs of the terms E2 f2 and E3 f3 

has the effect that the error influences caused by the lock-in 
amplifier are markedly reduced. The error. value given in 
equation (4.1.7) thus represents a high upper limit. 

Result: /69 

The error is less than or equal to 15% because of the 
measuring instruments and the measuring set-up. 

Measurements in small magnetic fields (300mT) lead to large 
mistakes . 

— The general adjustment specification can be given under -he 
condition E2 ( 3 < 8.66; 

B >/ l/4b (4.1.12) 


3 . 


1/ 0 L (e)~ U L l H) 

b 

’ U L < 3) 

3 ^ 


63 0 Vs/m. 1 


The third term 

r , ( 0 ) 

3 ,3 " 1 urn - u L rs) 


can be treated analogously to E 
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not attainable 
with measuring 
construction 


Figure 4.1.1: Adjustment specification for the magnetic 

field on the basis of acceptable errors and 
average mobility 

The 15% error boundary which can be reached with the present 
methods is too high to be satisfying. Physical errors would also 
have to be added to this value. For this reason steps must be 
taken which reduce error or at least offset errors. 

The procedure described in Chapter (2.4.3) lends itself to /70 
this goal. 
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Figure 4.1.2: Error offset by line approximation 

Mobility no longer results from a two-point measurement (S(0), 
S(B)), but rather from a successive measurement 

S ( 0 ) , S ( B ! ) , S(B 2 ), to S(B n ) 

4.1.3 Statistical procedure /71 

In Chapter (4.1.2) error for a single point in Figure (4.1.2) was 
estimated with the help of the concept of error sensitivity and 
with measuring instrument precision boundaries. These 
observations were meaningful and necessary. They indicated 
sensitive places in the measuring procedure and provided a 
preliminary adjustment specification for the magnetic field. 
Nevertheless, a universal error index is not possible since 
sensitivity factors are especially dependent not only upon the 
tests but also the working point. 

Already at the end of Chapter (4.1.2) a solution was suggested 
for reducing errors by line offsetting. The evaluation of the 
data made possible by such a balance becomes optimal with 
statistical methods: 

Input data control 

Statistical error balancing using a linear regression 
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Physical errors caused by parasitic resistances cannot be 
determined. If, however, the graph in accordance with Figure 
(4.1.2) deviates greatly from a straight line, this is indicated 
by large errors, which result from 

— The drift of the structural element 

Reaching the physical boundaries of the measuring procedure 

With that, the validity control of the evaluation equation 
recommended in Chapter (2.4.3) is carried out. 

4. 1.3.1 Control of input data /72 

The measuring value reacts to a change in setting dimensions by a 
transient process in the new stationary end value. 



r 


Figure 4.1.3: Transient process of a measuring value after 

changing the setting dimension 

Key: 1) Setting dimension 2) Measuring value 

Input data control should 

Increase the precision of input data by establishing median 
values from many measuring values 
— Eliminate measurement before a stationary end value has been 
attained, by controlling the spread of measuring values 
around the median value 

A storage bank is filled with k~measuring values of the same 
dimension. The measuring values are recorded at intervals T and 
deposited in the uppermost storage position of the bank. The 
lowermost (oldest) measuring value drops out. 
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Figure 4.1.4: Storage bank for k=4 

The following are formed from measuring values in the bank: 

— The median value 

-- The quadratic deviation from the median value, the variance: 

• k 

(4.1.13) 

1*1 

If the variance is too high, then the transient process is not 
complete and a new measuring value is taken into the pile. This 
procedure continues until the variance falls beneath a certain 
level chosen so that the spread around the median value is 
smaller than the error tolerance of the measuring instrument. 
Assuming normally distributed measuring values: /74 

* ~vt - 0 - * (4 - 1 ' 14) 

x : Median value 

u : Statistical threshold value 

a : Standard deviation 

U : True value 

The equation indicates in which interval around the median value 
the true value is to be found. According to German regulatory 
document DIN 55032 error should be estimated so that, with 95% 
confidence, the true value will not deviate more than 
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from the true value. For this situation /14 /: 


u = 1.96 

If we now set the following in the first 
k -*■ »■>) : 

f ‘*1 

a limit is reached with equations (4.1.14) and (4.1.15): 

r _ i <-U . i,36-6 (4.1.17) 

Tx ~ x x .{T 

Remarks : 

The approximation in equation (4.1.16) is relatively rough for 
sample test numbers k=4 - 6 used in the measuring device. Thus, 
it is not to be expected that error can be estimated accurately 
with equation (4.1.17). This procedure primarily provides 
datacontrol. Formation of the median value causes an error 
offset in measuring values and, for this reason, substantially 
improves the precision of the data input. 

4. 1.3.2 Error offsetting using a linear regression /75 

The input data prepared in accordance with the previous chapter 
allow us to make a diagram on the basis of Figure (4.1.2) The 
points have the coordinates: 



<*k : *k 1 


k - 1 , 2, . . . . , n 

A balancing line can be made using these points in a linear 
regression . 


(4.1.15) 

approximation (ideal for 
(.1. 1.16) 


y = ap + a^ • x (4.1. 18) 

A description of this procedure will be omitted here. Please 
refer to literature /14/. 

Results : 

Mobility can be found from the slope of the straight line: 
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Q 1 


b 


(4.1.19) 


The procedure produces two error criteria: 

1. Zero position error ag 

Ideally the balancing line goes through origin. If ag 
deviates substantially from 0, then we must assume drift 
attainment of the measuring procedure validity boundary. /76 

2. Coefficient error Aa^ 

With a statistical confidence of 95$!, the error determines 
the error interval around a^ where the true value of b lies. 

Both of these criteria determine the error index in the measuring 
room expression: 


Errors B/M 

X X 

The left-hand value is defined as the percent of intervals around 
the median value a^ = b: 



(4.1.20) 


The right-hand value means: 




(4.1.21) 


It indicates the zero point error in relation to the slope. This 
somewhat unusual definition became necessary because, unlike the 
foregoing error criterion, this one lacks a true reference 
dimension. Taking the absolute value alone is not permissible, 
since the dimension order of the zero point error is more 
difficult to recognize. 

4.1.4 Drift occarences /77 

The linear regression procedure can be optimized primarily by 
increasing magnetic field settings for every point of the 
mobility profile. The degree of freedom of approximation 
increases and error offsetting improves. This possibility does 
not fail due to pure time expenditure, but rather as a result of 
the drift in the measuring equipment and the test itself. 

The following considerations are aimed at eliminating drifts in 
the measuring room. The drift of the test cannot be influenced. 
It is indicated under the foregoing conditions in the error 
criteria according to equation (4.1.21). 
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The following table shows which drifts occur in the measuring 
device. The time period corresponds to seven magnetic field 
measurements per line approximation. 


Parameter 

0 

Drlftzeit 

s 

D r i 

2-) 1st 

0/ 

10 

f t 

3 ) Soil 
% 

U L 

245 

0,34 

0,1 

B 

35 

0,3 

1 

U GS 

245 

/ 

1 

u ns 

245 

/ 

1 

u gs 

245 

1,0 

0,1 


Table 4.1.5: Uncorrected drifts in the measuring 

construction 

Key: 1) Drift time 2) Actual % 3) Ideal % 

Actual values were attained from measurements in the measuring 
device. Ideal values are chosen so that drift does not go above 
the error boundary of the parameter error. The values and Ug a 
exceed permitted drift. It was found that the drift of 
(initial voltage of the lock-in amplifier) primarily is 
determined by the drift of Ug s . For this reason it is necessary 
only to make a correction in the alternating voltage of the gate. 


Transconductance is calculated from: 


5 = 


’ S g 

• u 3S • 40V 


(4.1.22) 


Before the correction it was assumed that the value of Ug s , set 
only one time, remains constant for all transconductance 
calculations during the measuring period. As such, the drift of 
Ug s had a direct and very sensitive influence on the measurement 
results . 


Correction: For every 

were recorded simultane 
exhibit a drift, howeve 
measuring . 


transconductance measurement Ug s and uj, 
ously. The voltage source continues to 
r, its influence is corrected by technical 


This precaution necessitates additional installation of a second 
digital multimeter compatible to the IEC-Bus, but it 7J3 

eliminates virtually the entire drift of the measuring device. 
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Parameter 

Drlftzeit 

s 

D r 
1st 
% 

f t 

Soil 

ot 

Jo 

U L 

245 

0,033 

0,1 

B 

35 

0,3 

1 

U GS 

245 

/ 

1 

U DS 

245 

/ 

1 

V 

0 

/ 



0,1 


Table 4.1.6: Corrected drifts of the measuring device 
• Key: 1) Drift time 2) Actual % 3) Ideal * 

The values in the table are related to seven magnetic field 
measurements. Since the ideal values are well satisfied, even a 
considerable increase in magnetic field settings on the side of 
the measuring device is not critical. The number of magnetic 
field settings can be governed completely by the type and quality 
of the test. In the case of tests with minimal mobility, they 
ought to be as large as the drift allows. 

4 . 2 Technical boundaries of the measuring procedure /8Q 

The measuring technique represented here has two limitations in 
practical use: 

Transconductance of the test decreases so much that the 
entrance sensitivity of the lock-in amplifier is exceded. 

Charge. carrier mobility b is so small that the magnetic 
transconductance change can no longer be evaluated. 

A numerical indication of the limiting values 

s S> b S 

essentially depends on which error boundary is allowed, 
evaluation involves an error indication for every value, 
operator may determine these limitations. Here we shall 
the causes and relationships of the limits and, finally, 
orientation values will be given. 

Regarding the t ransconductance boundary Sg : 

At the setting "max output stability", the lock-in amplifier has 
an initial sensitivity of 100 yV at full force. The smallest 
permissible value amounts to 20% of the full oscillation. 


Since 

the 

describe 
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(4.2.1) 


U L - 20 

From thie value the Sg for Ug a = 60 mV can he attained: 

S s = 400 (4.2.2) 

If the measuring device was put into operation approximately /81 
approximately one to two hours before the measurement, then the 
lock-in amplifier can be operated at the "normal" setting as 
well. As a result, initial sensitivity increases by a factor of 
100. At this low voltage, however, even the lock-in technique 
can no longer separate the noise portion with sufficient 
certainty so that the oscillation swings easily. A lower 
limitation which can still be evaluated by technical measuring is 
indicated with: 


Sg = 20 jjS 

Regarding the mobility boundary bg 
Equation (4.1.4) exhibits an additional term 


. r A J_ D, (o) . .. 'iVu 


(4.2.3) 


(4.2.4) 


which describes the error portion of the 
help of the equation 


b 


_L. 

B 



lock-in 


amplifier. With 


(4.2.5) 


a limit bg can be indicated. The error portion of the lock-'in is 
set at 20*. The measuring instrument tolerance amounts to: 


MJl 

Ul 


0,01 


With these conditions the following is attained 
(4.2.4) : 

U, lo) 

U,(0)-V t (3) 



for equation 


(4.2.6) 


The initial lock-in amplifier voltage decreases by 2.5* in the 
magnetic field. This small change is the useful information! 

/82 

u L (0) = 1.025 u l (B) 
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With this the following is approximately true for equation 
(4.2.5) : 


b = b s - 


For a strong magnetic field (B 
expressed by: 

b s = 

Measurements on test III-49 (see Figure (5.1.1) to Figure 
(5.1.3)) verify this value. The graphic representation of the 
profile shows, however, that even below this value a mobility 
profile with a steady course can be obtained. Error offsetting 
of the linear regression assures that error in the end result 
lies below occassionally high parameter errors. Since the values 
spread widely, the necessary statistical certainty of 95* can be 
maintained only in a relatively large error interval. 

These observations can be used to indicate the technical 
boundaries : 

S >, 0.02 mS 
b 1760 cm 2 /Vs 

4 . 3 Setting recommendations for parameter measurements / 1 * 

The parameter settings 

u gs* U DS ' B 

essentially determine measuring precision of the end result. 

They should be set so that the physical errors and the technical 
measuring errors maintain the same dimensions. This brings about 
end results with minimum errors and little effort. Occasionally, 
both of these error sources cannot be reduced simultaneously, so 
that the setting recommendations provide a suitable compromise 
between : 

1. Ug s minimum 

2. Up S minimum 

3. B equally distributed 

4. B maximum 

5. Up lock-in at full amplitude 



= 900 mT) the lower limit of b is 
1760 cm 2 / Vs (4.2.8) 


♦Section 4.3 is included at the very end of the original text, 
with page numbers as indicated herein. 
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The compromise between 3. and 4. is found with equation (4.1.11). 
There the lower limit of B is: ^ 

Vb 

b : Estimated value of mobility 

This condition is quite astutely described. In practice the 
following demand proves to be sufficient: 

1 (4.3.1) 

Bs Tb. 

The estimated value of b can be taken from expected test data: 


Test 

b / cm 2 / Vs 

Bmin/ 1 ^ 

Epitaxy 

3000 

480 

Diffusion 

2500 

580 

Ion implantation 

4000 

360 

HETERO 

6000 

240 


Table 4.3.1: Estimated value of mobility in various tests 


The remaining value range of B ra ^ n up to 999 mT is carried out in 
similar fasion. However, it is not practical to set B-values 
under 500 mT. 


The compromise between 1., 2. and 5. is found when the values of 
Ug s and Ups bring about full amplitude in the lowest range of the 
lock-in. The following table shows the relationship between 
transconductance test dimensions as well as Sg (entrance 
sensitivity of the lock-in amplifier) and Ug s . 


S = 


4 ^ 

' loV 


f(u gs* S E ) /mS 


(4.3.2) 
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U QS 

mV 

Se 

10 yV 

nsitivtty 
30 yV 

100 yV 

30 

0.06 - 0.39 

0.18 - 1.18 

0.61 - 3.94 

40 

0.05 - 0.30 

0.14 - 0.89 

0.45 - 2.95 

50 

0.04 - 0.24 

0.11 - 0.71 

0.36 - 2.36 

60 

0.03 - 0.20 

0.09' - 0.54 

0.30 - 1.97 

70 

0.03 - 0.17 

0.08 - 0.51 

0.26 - 1.69 

80 

0.02 - 0.15 

0.07 - 0.44 

0.23 - 1.48 

90 

0.02 - 0.13 

0.06 - 0.39 

0.20 - 1.31 



S„ - 





u 0 



Table 4.3.2: Calculable transconductance range with 

dependency of Ug B and Sg 


Transconductance boundaries are defined so that! 


Z2 


S = S u : U L = 2 V 
S = S 0 : U L = 13 V 


(4.3.3) 


U L ! Exit direct voltage of the lock-in amplifier 

With program "B-transconductance measurements" we establish a 
transconductance profile 


S = S(U QS , U DS ) 

From this profile we seek the measuring series S(Uq s , U ds = 
const.), which at the smallest possible Upg and Ug s fits into 
value range 'S u - S 0 ’ of Table (4,3.2). 

For partial measurements 


— Ug a at the gate contact (Uq S > 0 V) 

— UdS on the substrate (Ugg U^) 


must be chosen to be very small. Otherwise both voltages should 
be of the same amplitude (standard value). 
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5 


Results 


783 


This chapter contains a selection of measuring results made with 
the MAGNETO-TRANSCONDUCTANCE-MOBILITY-PROFILING procedure. 

MESFETS of various technologies were investigated with a wide 
value range for charge carrier mobility. 

The mobility profiles represented are complemented by remarks on 
the technology of structural element and the measuring 

parameters. The selection of measuring results covers the 
characteristic strengths and weaknesses of the measuring 
procedure . 

5.1 Diffused MESFETs 

The active layer is produced by diffusion in semi-isolating QaAe. 
Diffusion time and temperature are decisive technological data 
for the concentration of the charge carriers. A look at the 
comparison between measured results and technology of two 
diffused MESFETs which were investigated clarifies the influence 
of technology on the structural element. 

a) III-49 

The measuring results are represented in Figures (5.1.1) to 
(5.1.3) 

Technical measuring remarks: 784 

Charge carrier mobility lies in the lower boundary range at 
approximately 1700 cm^/Vs. Nevertheless, a steady mobility 
profile was obtained. The b-value of U GS = -1.7 V falls out of 
the path. In Figure (5,1,2) the triangles represent magnetic 
field settings for this value and illustrate that apparently 
drift occurrences lead to this measuring error. The error course 
of this profile (see Figure 5.1.3), like all the other profiles, 
agrees well with the expected values. The very high values in 
the proximity of the boundary layers indicate the boundaries of 
the measuring procedure. They have their basis in high drift 
occurrences (boundary layer effects) and lower mobility values. 

b) III-78 

The measuring results are represented in Figure (5.1.4). 

This profile represents a portion of the entire range. Both 
profiles will be compared following the figures. 
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Mobility Profile May 16. 1983 /86 


Teat 

III-49 

10.7 

Data input 

BIII49 

10.7 T 

Measuring time 

273 minutes 

Sensitivity 

100 UV 


Magnetic field /MT 

600 / 

970 / 90 

Alternating voltage 

50MV / 

1KHZ 

U ( DS ) 

.1 V 




U<GS> 1 

BE 

1 ERROR 

1 
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1 


V 1 
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i y. 

1 

y. 

1 
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1 
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1 

I 


+0. 40 1 
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1 


+ 0. 00 1 
1 
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1 

— , ♦*« 
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-0.10 i 

1 860 
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1 5.9 
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-0.20 1 
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1 4.7 
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- .9 
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-0. 30 i 
1 

1220 

l i o 
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1 

• 4 

! 
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1 3.1 

l 

1 

. 5 

1 

\ 

i 
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1 
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1 

1 900 
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1 

1 
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I 
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1 
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1 

1 

- . i 

1 
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-0.80 1 

1900 

1 1.5 

1 

- .i 

1 


• -0.98 1 

1 

1 850 

1 1.8 

1 

1 

* — 

1 

1 


1 

-1.08 1 

1830 

1 1.5 

r 

i 

- .0 

1 

1 


-1.10 I 

1828 ■ 

1 2.3 

i 

+ .5 

1 


-1.20 1 
1 

1 7S0 

I 2.4 

I 

i 

i 

- .4 

1 

1 

1 


- 1 . 38 1 

1770 

1 .5 

l 

f 

“ . 0 


-1.40 1 

1740 

1 3.3 

i 
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1 


-1 . 50 ! 

i 

1688 

1 1.9 
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i 

- .2 

1 

1 


-1.60 i 

1640 

1 8.2 

I 

i 

- 1.6 

1 

1 


-1.70 1 

1820 

1 54.9 

i 

- 1.4 

1 


-1.00 f 
| 

1500 

1 23.3 

1 

i 

1 

- 4.8 

1 

| 


-1.90 1 

1460 . 

1 14.5 

l 

i 

- 3.0 

1 

1 


-2. 00 1 

1340 

1 51.9 

i 

-10.7 

1 


-2.10 1 
1 

1320 

1 6.9 

1 

i 

i 

- 1.2 

l 

1 


-2.20 1 
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l 

i 

-14. 1 

1 
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B'igure 5.1.3: III-49 Tabulated mobility profile 
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3000 f 


2000 f 
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.5 

Figure 5.1.4: 


u ( 9 s ) = 50 mV 

U ( DS ) = 100 mV 


0 -.6 - 1.2 - 1.0 

U(GS)/V 

III-78 Mobility profile 



1 1 1-49 

1 1 1-78 

N D /10 18 cm' 3 

T Diff 

t Diff 

A) Gateabsenkung 

W c, " 2/Vs 

^max 

1 ,5. . .2.5 
800 °C 
10 min 
B) gering 
1800 
15 GHz 

0,6. . .0,7 
900 °C 
36 min 
c) grod 
4000 
21 GHz 


Table 5.1.5.* Data comparison for two diffused tests 

Key: A) Gate decrease B) Small C) Large 


The different concentrations can be attributed to diffusion data 
in connection with the gate decrease. The decline in 188 
mobility at the gate contact has different causes in the tests: 






III-49: Surface effects 

III-78 : Steady increase of the disturbance concentration 

The data of III-78 are unusual for a diffused MESFET and 
particularly in relation to its high mobility "'nd high boundary 
frequency very interesting. According to th>' 4 structions of the 
producer of the structural element, the progress of the charge 
carrier concentration can be closed upon depletion of the doping 
material source during layer production. 

5 . 2 Epi t axial -MBS FETs 

a) I 208-b 

This structural element was produced with liquid phase epitaxy 
(LPE). The measuring results are represented in Figures (5.2.1) 
and (5.2.2). 

Technical measuring remarks: 

In spite of the small number of four B-values per line 
approximation, a very precise mobility profile could be obtained 
(Figure 5.2.1). A comparison of the error criterion "B/M" with 
test III-49 clearly shows that the diffused MESFET exhibits a 
higher drift and that, at nearly identical test mobility, this 
drift is responsible for the poorer results of the diffused 
tests. In observing the absolute values of the error criterion 
"errors in percentage" it should be remembered, that the 
statistic errors have been evaluated there. A physical error of 
5% maximum would have to be added in. 

The error boundary increases sharply at I 208-b in the /89 

proximity of SI. This allows us to conclude that there are 
technological difficulties in the boundary surface. 

b) ALF 1003/1 

This MESFET, produced in gas phase epitaxy (VPE), was 
investigated as an example for an industrial model. 

In Chapter (5.1) both of the diffused tests 111-49 and III-78 
were introduced. At absolute value and course the profile of the 
I 208-b has e mobility profile similar to that of the III-49. 

This is also true for tests ALF 1003/1 and III-78. In the tests 
mentioned last, the decline of mobility at the gate contact can 
be attributed to the steady increase of charge carrier 
concentration in both cases. The other tests each exhibit high 
concentrations, which lead to a decrease of charge carrier 
mobility. 
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MOBILITY profile 


APRIL 28, 1983 


TEST 

DATA INPUT 
MEASURING TIME 
SENSITIVITY 
MAGNETIC FIELD /MT 
ALTERNATING VOLTAGE 
U ( DS ) 


1-208 B 0.5 
BI208B 28.03 T 
85 MINUTES 
30 UV 

601 / 900 / 10G 
50MV / 1KHZ 
.05 V 


LKCiS) 
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Figure 5.2.1: I 208-b Tabular mobility profile 
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Figure 5.2.3: 









5.3 


Hetero-structure MESFET 


/92 


The construction of a hetero-MESFET differs considerably from the 
homogenously doped structural elements considered thus far. The 
following figure shows a cross sectional view of the structural 
element : 


Electronic 
gas ' 


n • Aly Co^.y Aft 


und. Aly Go yy AS 


GoAs 

( undoped ,'p’l 


s. I. GoAs substrate 



Figure 5.3.1: Structure of a hetero-MESFET 


The active layer was produced using molecular ray epitaxy. The 
literature /17,18,19/ discusses the functioning process of the 
structural element. The most important characteristic is the 
formation of a "two-dimensional electronic gas" in the upper 
boundary layer of the undoped (p“) GaAs. This produces a very 
high mobility at simultaneously high charge carrier concentration 
in a very thin layer. The basis for the high mobility is the 
very low concentration of ionized disturbance place atoms and the 
homogenous grid structrue which results from .it. 


The mobility profile could be determined only insufficiently 
using integral measuring methods. Thus the Hall measurements 
yielded a charge carrier mobility of 4500 cm^/Vs. With the /93 
following measurement the actual mobility and its course are 
recorded for the first time and with it a technical measuring 
substantiation of the "electronic gas": 
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Figure 5.3.2: 


a • 24.apr 83 
■ : 1 1 . oa| 83 

D-6-9a Mobility profile 


U ( GS 3 /V 


The research on this structure has not yet been developed very 
far. This measuring procedure may contribute a great deal to the 
further development and understanding of the structural element. 

Remarks on the measuring technique /94 

In Figure (5.3.1) the profiles of two measurements have been 
superimposed. The high degree of agreement proves that the 
measuring procedure can be reproduced successfully. Median 
deviance between the measurements amounts to 2.3%. Maximum ' 
deviance (at SI) amounts to 9%. 


6 Summary /95 


This measuring procedure enables us to make an actual 
determination of local mobility in the channel. This procedure 
is applicable to all current technologies of short-channel- 
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MESFETs and provides up to b^slOOO cm^/Vs useful data. 
Technological inferences regarding surface disturbances, grid 
voltage, and quality of the boundary layers are made possible by 
means of comparison with other parameters in the channel. 

Measuring parameters Uq s and Ug 8 make it possible to set and 
cancel depth and layer thickness. This is not possible with any 
other procedure. As a result of this study, profiles of high 
precision over the entirety of the channel have been made 
available. The measuring results for both course and absolute 
value stand up against results of other works (/1,4,5/). 
Furthermore, this study provides the first mobility profile of 
modulation doped hetero-structures measured using this procedure. 
The results provide a very good look at the functioning processes 
and quality of this new type of structure. 

Other focal points of this study lie in the areas of 

Measuring room automation. 

— Error control and prevention 


Measuring room automation has grown with lEC-Bus compatibility of 
measuring equipment and instruments. Further progress in this 
direction is possible only if a lock-in amplifier is developed 
whose initial sensitivity can be controlled directly by the 
measuring program. The measuring room can be serviced by a 
trained lab assitant. /96 


The author knows of no publications which deal with error 
observations of the measuring values. The error criteria 
developed in this study provide results which agree very closely 
with theoretically expected procedures. From the observations of 
statistical and physical errors, a value of 

f < 10 % 


can be given as 
cwr/ Vs . 


an error orientation 


for mobilities over 3000 
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b-roeasuremen t 


11.05. 1903 


Page: 1 




100 rem ### 

110 rem ### be B ftiuigsla ### 

120 rem N«N 

130 open4,4sopenl6, 16s ch* c chr» < 147) sprintch#: print sprint 

140 print #16, "b8"ssd* s ‘"s3v"sc*""c99"i l» B "199"sr**"xr" lprint#4, ’'s3v0000xr 11 

150 ;<z»292.2sxr m 430.3sxb B 89:rem threshold value for s®95’/. 

160 print" data inputs "sprintiprint 
170 print: input" dote ";da* 

100 printsinput" t««t "jpn* 

190 printsprint" u(qs) "sprint" " 

200 input" alerting value v "jga: input" end vollJe v "jge 

210 input" atep width v " j gs : gf ■ i nt < <ge-ga ) /gs ) +2; gn=0 

220 ge=ge+.001 Igs: i f g f < 1 then gosub 1030: got cl 90 

230 it <ga-. 7) * (ge-. 7) <0 or (ga+5. 5) * <ge+5. 5) <0 then gosub 1090igotol90 

240 print sprint" b "sprint" " 

250 input" « tartlng value mt " > baS input" end value mt " j be 

260 input" atep width m t "jbs:ba=ba-.01*bs:bf«int < (be-ba) /bs>+3 

270 if bf<l then gosub 1030sgoto240 

200 if gftbf>300 then gosub 1060:gotol90 

290 if be>1000 or b a. '>1000 then gosub 1090s goto240 

300 printsinput" u<ds) v ";ud 

310 if ud>.3 then gosub 1090sgoto 300 

320 printsprint" aitcmeting voltage "sprint" — " 

330 input" frequency kh: ";f 

340 input” amplitude mv a: af *=str» (a) +"mv / "+str* <f ) + "khz * 

350 printsprint" aeaeuring range "•print" " 

360 input" sensitivity uv "jse 

370 input" dc-vol tmeter v ")br:br“20000/br 

300 input" ac-vol tmeter mv "iarsar != 20000/ar 

390 printch*sprintsprint:print" aagnetic field nomoiiiation "-gosub 1120 
400 dim s (gf , bf ), a*( 152) : ww=0t ti $="000000": print " ■easureaent " 

410 for gaga to ge step gs: g=int (10000*g+. 5) /10000s gosub 620s rem r?<s apga .'ourco 
420 gn B gn+lsbn=Osprintsprint" u (gs)'*" | gj " v"sprintsprint 
430 for b=be+bs to ba step -bs 
440 bn=bn+l 

450 if bn>l then gosub 1300sprint" ";b,s rem b-adjuetaent ' ( 

460 if bn=l then print" b /mt",,"s /ms":printsprint" ~6", 

470 gosub 030s rem k-dmm 10+24 

400 if bn>l or sgn < (13-r) t (r-2) >=l then 530 

490 printsprintsprint" korrektur der lock-in eingangsempf ind- lichkeit" 

500 printsprint: input" sensitivity uv ";se 

510 printsprint" u<gs)=";g; " v"sprint:printiprint" "jO, 

520 goto 470 
530 next b 

540 print#4, "s3v0000xr":wz=240: gosub 790:if ww=0 then ww=l:bn=0:goto430 
550 print#4, "s3v0000xr"sprint:print:ne>:t g 
560 m$=ti tsprint#16, "b8" 

570 print#4, "s3v0000xr" sprint sprint"- end of .aeaaurcBent/ evaluation" 

580 printsprintsprint": .atore delate ("sinputme* . 

590 if me$="s" then gosub 1860: goto610i rem 3000 ,toro 
600 if me*<> "1" then 580 . 

610 close6;c]ose4;closel6:end 
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b-measuremen t 


11. '05. 1903 


Page: 2‘ 


620 rem IH»*» rS-.s voitagn «ourc« MOM 
630 open 18, 9 

640 ua* a "ual "+str* <g) I ub#**"ubl "+str* (ud> :printM18,ua*sprintM18,ub* 

650 wz a 60sgotub 790 

660 inputMlB, au*s an«l en (au») s i f an*24 then an**an-l 

670 ua * val(right*(au$,an-3) ) sub ■ val (right#<au*, an-15) ) 

680 if uaOg or ub <>ud then pri nt .. voltaiB , 0U rce an" sgosub750sgoto640 
690 printM16, "aB" 

700 cl osel8 jreturn 

710 rent MOM input loop ■ 1 MO# 

720 get j#sget'j*» get j* s get j*i j*=" " s print : pri nt" 'j/n'?" 

730 get j«i if j«-“" then 730 

740 print, jMsreturn 

750 rent MOM input loop 2 MM# 

760 get j*sgetj$s get j»i get j#s j*«"" spri nts print" 'j'?" 

770 get j#s if then 770 

780 print, j*»return 

790 wt=»ti srem MMM waiting loop MOM 
800 if wt+wz > ti then 800 
810 return 
820 rem MMM 

830 rem MMMmeasuring value romkei thl ey multimeter 
840 open7,24:open20, 10i j=0swz a 300sgosub 790 
850 inputM20,w*J if stOOthen 850 
860 inputM7,z*s if stOOthen 860 

870 j=j+ls ifw* a "over"then printsprint"overflow"sstop 
880 if z* a "over "then printsprint"overf low" j stop 
090 r5=r4sr4=r3ir3“r2sr2=rlsrl=val <w*)/br 

900 z5 a z4s z4 a z3s z3 a z2sz2 a z 1 s z l a val (z*) /ars wz»150s gosub 790tif j<6 then 850 
910 if j > 1 5 then printiprint" Idool value dan not convergent " J j a 8sstop 
920 r a <rl+r2+r3+r4+r5) /5s rq« <r l-r > 62+<r2-r > 62+<r3-r> 62+<r4-r )62+ <r5-r) 6 
930 z = <z 1 +z2+z3+z4+z5) /Si zq«= <r 1 — z ) 62+ <z2-z ) 62+ <z3-z > 62+ <z4-z ) 62+ (z5-z ) 6 
940 if zq a 0 then 960 

950 if z/sqr <zq/4X>!2 then 850s rem xz s 957.wk. for f< 0.37, 

960 if rq=0 then 990 
970 ra=sqr <rq/4> sr>i a r/ra 

980 if r>:<xr then 850i rem xrx 95%wk. for f< 0.27. statement from lopiace 
990 s<gn,bn)= r /z tse/ 1 Oi rz=int ( 1000<s <gn, bn ) +. 5) /1000 
1000 print, rz 

1010 if sign, bn) >s<gn, 1) then print, "drift error" i stopib=0: gn=gn-l j g a g-gs 

1020 clO5e7:close20ireturn 

1030 rem MMM sr error input 1 MMM 

1040 printiprint" consult instructions "sprints wz=60sgosub 790 
1050 return 

1060 rem MMM sr error input 2 MMM 

1070 printiprint" amount of too high "sprintswz»60igcsub 790 

1080 return 

1090 rem MMM sr error inputs 3 MMM 

1100 printiprint" measuring range exceeded " s pr i n t : wz=60s gosub 790 

1110 return 
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b-meaaurement 


11.09. 1963 


Page : 


3 


1120 rem HHH ■■■■ur!n< ran»s, b, d *f(b) ### 

1 130 printsprint" nax'au* vaiu** for b*b <max) "t print 
1140 mb*1000s if be<95 then mb *100 

1150 printlprintjprintjprint" BeBaurinir ran<«> tola Uttar opan 11 

1160 print sprint sprint sprint" b <ma>: ) ■" jmbj " mtesla "iqosub 750 

1170 d* B "0Q0Q" s pr i nt*4, sd*+d*+c*+J *+r* 

1160 wz*60sgosub 790 

11.90 print#16, "b7"sgo»ub 790iprint#16, M b7"igoBub 790 

1200 if d**"0000" then print#16, "a6":gasub 790:print#J6, "a6"igosub 790 
1210 print#4,sd*+ M 9?99"+c«+l*+r* 

1220 printch*sprint" talker"sprirttsprint" taste inhibit on ?" sgosub 750 
1230 printch** print: print sprint" b* "|mbj" mtesla over voltage iourc." 

1240 print}print:print»print" b*"(mbj" mtesla ??!"jgosub 750 
1250 printch*: printiprint" taste inhibit out ?" igosub 750 
1260 print sprint#4 , sd*+"00Q0"+c*+l *+r* 

1270 a 1*8. 6063s a2*l . 3937e-3j def f n us <:< )*i nt <ailx+a2*>;o2+. 5) 

1280 deffn ur<yl=al+2ta2*y 
1290 return 

1300 rem *## b-odjuatnent 

1310 d0=fn us(b> 

1320 cpenl2, 12idl*d0 
1330 z=0 

1340 d*=r ight* ( "0000"+str* <d 1 ) * 4 > 

1350 aa**sd*+d*+c*+l *+r* 

1360 J=0s print#4, aa* 

1370 wz = < 1-sgn (z > ) *240+120: gosub 790 
1380 input#12,b*s if st <> 0 then 1380 
1390 if 1 en <b*) <6 then 1330 
1400 wz*10:gosub 790 
1410 er*=right*(b*,4) 

1420 e«val <er«> /10*mb/ 1000 
1430 el*=l ef t* (b*, len <b*>-4) 

1440 if el*="ee " or el*="6e “ then e=e sgoto 1470 

1450 if elf^'eee" or el* 3 "6ee* then e=-e sgoto 1470 

1460 if el*="e " or el5="6 " then e=e+1000 sgoto 1470 

1470 if el**"e " or el*="6 " then print, , "overload"sgoto 750 

1480 e4=e3:e3=e2se2=el:el=es j*j+l: if j<!5 then 1380 

1490 bb=(el+e2+e3+e4) /4: bq=<el-bb)62+(e2-bb>o2+(e3-bb) o2+<e4-bb)o2 

1500 ab=sqr <bq/3) sbx«bb/ab 

1510 if abs<bbK.02*mb then 1130 

1520 if b>!< xb then 1380: rem 957. wkt for f<l7. (atatasent from ia P ioc«). 

1530 rf=abs( <b-bb) /b! : i f rf'1.01 and z>0then 1590 

1540 2=2 + 1 S if2 >10thenZ*5!prints print" ‘deni value unattainable "sstop 
1550 dd=int (dl- (bb-b)lfnur < <b+bb)/2) +.5) 

1560 if dd>9999 then dd=9999 
1570 if dcj< 1 then dd=5 
1560 dl=dd:goto 1340 
1590 closel2sreturn 
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1600 rem ### ■‘■'•r* ### 

1610 printch*:print:print") lo . d t«bi« of content* "iqoto 1840 

1620 printiprint "opening of document" * 

1630 printjprint"/ !'•»•« ■or# than 30 block* free7“ *'jgOSltb 710 
1640 if J**"n" then print -now diskette" stop 

1650 print" from teet number end knowledge Input new document nemo" 

1660 print:print"namo» <(15)" 

1670 printi input" <1 o , : u "« t ' t "■■•"jdn* 

1680 if len(dn*)>15 then prinf’^ocument n««c longer then 1G cher«Ctere"'!goto 1650 
1690 n*= M 0"+ln*+"j "+dn*t-",«,H l ’j open 100,8, 10, n* 

1700 f 2«14» i f bf >14. then f2»bf 
1710 dim s* < f 2) 

1720 s* (1 )=da*{ s* (2)*pn*i s*<3) a str*(ga) i 5* (4 ) a str* (ge) is* (5) =str*(gs) 

1730 E*(6)“str*(gf ) i c* (7) -str* (ba*. 01 tbs) j s* (8>“str*(be) i s* (9) *str* <bs) 

1740 s*(10) a str*(bf > i s* ( 1 1 ) a str* <ud) i s* < 12) «=af *: s* < 13) a str* (se) i**(14) a mt 

1750 for m“l to 14 

1760 print# 100, s* (m) ,* chr* ( 13? J 

1770 next 

1780 for n a l to gf 
1790 for m a l to bf. 

1800 s* (in) a str*(s (n,m) ) i print #100, sf (m) iehf*( 1 3) j 

IB 10 next m 

1820 next n 

1830 close lOOtreturn 

1840 rem ### loed table of content* ### 

1850 k«0idim a«<152) 

1860 printiprintspr int i printiprint 
1870 input"laufwerk 0/1 "jin* 

1880 printch*;print:printsprint" pienet wait 30 *econd* " 

1890 openlS, 8, 15 

1900 print# 15, "m-e"chr* (212) chr* (237) 

1910 open2, 8, 0, M *"+ln*+"j "+"b»" 

1920 get#2, a*: get#2, a* 

1930 1»1 

1940 get#2, a*: get #2, a* 

1950 a=0 

1960 get#2,a*:if a*<>"" then a a asc(a*> 

1970 get#2,a*:if a*<>"" then a »a+256*asc (a*) 

1980 a* <1 ) =str* (a) + "" 

1990 get#2,a*ia*(l >*=a* < 1 ) +a*» if a*<>" " and 1<>1 then k»l 
2000 get#2, a*; a* (1 ) =a* (1 ) +a* 

2010 if a* <> then 2000 

2020 if kOl then l=l+l» goto 1940 
2030 cl06e2:closel5 

2040 rem ### monitor Vfcri,. ,n of teblo of content# . h## 

2050 print ch* 

2060 for n c 0 to 1 
2070 printa*(n) 

2080 next 
2030 goto 1620 
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100 rum MN 

110 rem ♦!## voraion lie-profile HH4 

120 rem H#»i 

130 dim »l»{14)ich* ,l! chr«(147) 

140 printch*! print* print" road n uoaaureeent P ragree»icn tram diskette" 

130 print l printiprint "choice of »OM»urtng progreeelbn" 

140 pr i nt S pr i nt " collect table of contenta of dlakette t " I goBLib 1010 

170 printiprintiprintipnntipnnt 
180 input" 2 document n»»«")dn* 

190 n*»ln«*"i "♦dn**", s,r"»opnn 90,B,9,n* 

200 printch*sprintiprint" »i 0 ed protocol bead" 

210 for m B l to 14 
220 inputH90, si* (m) 

230 next 

240 f l«0l prints'. " *obllUy profile ^ "(filSU) 

250 print" printiprint 

240 print" teat H |Slt(2) iprintiprint" alternating voltage "j»l*(12> 

270 printiprint" be»f(u(gs>) /cm2v-ls-l":printiprint 
2B0 print, "u(ds> * "jsl*(ll)|" Vtprint 

290 printiprint" if ■•■■uring progreaalon correct "jgosub 930 j 

300 if J* *"J" then 330 . • 1 

310 close90iprinti print" now neaauring progression" ; gosub 930t i f j**" j" then 140 

320 goto 920 s 

330 gf **val (sl*<6) ) ibf »val (si* < 10) ) 1 4 l»gf : i * bf>gf then fl»bf 

340 dim x(fl)idim y(fl)trem for plotter 

350 dim s(gf ,bf > , b (gf , 3) 1 1*395/ <bf — 1 ) + 179: rem o PP rox. threehoid voiue s*957. bs 
340 for gn «*1 to gf 

370 for bn »1 :o bf f 

380 input#90,Ee*iB<qn,bn)' E val (ea*) iprints (gn,bn) j 

390 next bn ! 

400 next gn ‘ 

410 close90 

420 for gn *: gf • 

430 if sign,: s :hen 490 

440 for bn »L ' o4 j> 

450 if sign,; then 470 r 

440 sign, bn! =i . (gn, 1 ) /s (gn, bn)-l ) I 

470 next bn 
480 next gn 

490 ga a val(sl» tge»val (si* <4) ) igs^val (si* (5) ) i gn-0 

500 ba»val<6ll tbe^val <sl*<8) ) :bs=val < s 1 * < 9 ) > ibn*0 

510 for g=ga t ♦ step gs 

520 gn»gn+lixc s a 0i xs«0i xy=0i bn 3 l 

530 for b=be t b . step -bs 

540 bn=bn+lixc c- (b/1000) 62:x5 3 XB+b/1000: ys 3 ys+s <gn, bn) i xy 3 xy+b/1000ls (gn, bn) 

550 next b 

540 nn=<bf-l)» :-.« 602 ib (gn, l) a < <bf-l) *xy-xs*ys) /nmb (gn,2)=(xq*ys-xsfxy) /nn:bn=; 

570 sq=b(gn,2>62 

580 for b-be to ba step -bs 

590 bn=bn+l i cq=sq+ ( b ( gn, bn) -b (gn, 2>-b (gn, 1 > tb/ 1000) o2 
400 next b 

410 sa~sqr (sq/ (bf -3) ) 

420 b <gn,3)’=‘ttsa*sqr < <bf-l ) /nn)/b(gn, 1) :b (gn, 2> a 100*b(gn,2) /b(gn, 1) 

630 next g 

640 printiprint, " . Ubl °" ,, "pi otter "igosub 950 I 

650 'if j* 3 "p" then gosub 1310igoto920 
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660 open6, 6s f l»0sprint#6, , " 'mobility promo "|#1*<1) 

670 print#6 f , " »n»«.«eMs«=mB«Bmi»«i.«Ba=MBs>«««M««*BB«*BB>*jp r i r ,t#6sprint#6 
680 print#6,," test "jsl*(i .print#6, f document name "jdn* 

6?0 z**»sl*<14)si*»~str*<val (lef t* <z*, 2> ) «60+va) (mid* (z*, 3, 2) ) ) + " min" 

700 print#6, , "aeosurintf time "jZ* 

710 pri nt#6ipri nt#6, , "sensi ti vi ty "jb1*(13)j" uv" 

720 print#6,, "magnetic field /mt "j int <ba+.01 tbs> j "/"jbej jbs 
730 print#6sprint#6, "alternating voltage "|sl*(12) 

740 print#6 f> , "u(ds> ‘ si* < 1 1 ) } " v"s pri nt#6i pr int#6 

750 print#6sprint#6sprint,#6 f , "ii u(gs) X be X «rror X b/m S" 
760 pr i nt#6, , "<a! v X cmo2/v3 ini 7. £ V. £" 

770 q2*= "ACCCCCCCCCCaCCCCCCCCCCCC£CCCCCCCCCaCCCCCr.cna" sprint #6, f q2* 

780 f 1*=" a a a a a" 

790 1 2t=" a 52.79 a 99999 a 999.9 a s99.9 a" 

800 clos56sopen2,6,2iopenl,6, 1 s c*="£"+chr* (29) 

810 gn=0sprint#2, f 2* 

320 for g=ga to ge step gs 
830 gn=gn+l s f l=f 1 + 1 

840 pr l nt#l ,c*;g;c*; int ( lOOOtb (gn, 1 > + .5) *10j c* jb (gn, 3) +.05J c*}b (gn, 2) +, 05 1 c* 
850 if fls'3 then 870 

860 f 1 =0: print #2, f i*» print#! , c*;c*}c*{c*jc*sprint#2 f f2* 

670 next g 
E80 close6 

090 printchtsprintiprint" plotter expression 11 : gosub 930 
700 if j*="j" then gosub 1300sgoto 920 
910 if j*<> "n" then 890 

920 printch*:print:print: pr int" enrie"send 

930 rem input loop v 1 

940 printsprintsget j*sge< jfsget j*s j*="" sprint" 'j/n'" 

950 get j*i if j*="" then 950 
960 return 

970 rem input loop 2 

980 printsprinnyeijmsget j*sget j*i j*="" sprint" ’j" 

990 get j*s if j*0"j" then 990 
1000 return 

1010 rem ft## loud table of contents * ### 

1020 if yyOl then 1030 

1030 printsprint" new starting letter * "igosuh 930sif j*="j" then 1070 

1040 if yy=l then 1240 

1050 k=0:di ma* < 152) s yy=l 

1060 printsprintsprintsprintsprint 

1070 input" drive 0/1 "flntlprint: input" input "startihg letter of document" .*b*. 
1080 printchtsprintsprintsprint" picose wait 30 seconds " 

1090 openl5 t 8, 15 

1100 print#15, "m-e"chr* <2121 chr* (237) 

1110 open2,B,0, "*"+in*+"s "+ab*+"t" 

’1120 qet#2 f a*sget#2,a* 

1130 1=1 

1140 get #2, a*s get.#2, a* 

1150 a=0 

1160 get#2,s»sif a*<>"" then a =asc(a*) 

1170 get#2,a*sif a»<>"" then a =a+256* asc (a*) . 

1180 a* <i ) =str* (a) +" " 

1190 get #2, a*sa*(l)=a*(l) +a*s i f a*<>" " and K>1 then k=l 
1200 get #2, a*sa*(i)=a*(l)+a* 

1210 if a* <> "" then 1200 

1220 if k<>l then l=l+ls goto 1140 

1230 k=0sclose2sclosel5 

1240 rem table of cont. on monitor 

1250 printch* 1 • 8 4 

1260 for n*0 to 1 
1270 prints* !n) 
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1300 rem ar #K# digi plot ### 

1310 printjprintjprint" function : "iprint" CCCCCCCC'i print 
1320 input"be B f (u (gs) ) - aqr ( )«f(b) - u6b "jdw* 

1330 if dwIO'u” and dw«<>"b" and dw^O**" then 1320 
1340 xl*400iyl B 200jy2 B 1500irem diugro» origin coordinate* , 

13S0 nM B 4my B 5ldX B 225ldy B 180l rem length and nu»ber of axi* mcUom 

1360 ad B 150« sg B 3« sp B ag+l ihb B B0*7i rem sdtabst. d-s / sg: s.gr. /hbs 5-brei te 

1370 pr.>ntch*Jprintiprint" aaxiaua value of be" 

1380 input" cm2/vs "jum:um B int < <um-l) /1000+11 t 1 000 ibm B um/ 10000 

1390 if dw* B "b" then 1460 

1400 printchtjprintiprint" •cole for^ u(gs)":print 

1410 print", 5 0 -.5 ...-1.5 v 0,5":print" 1 0 -1 ...-3 v 1.0" 

1420 input" ";uw 

1430 printjprintiprint" nxi* acaie for "jprint 

1440 print" be(max) ="|umj "cm2/vs":print" be(min) B 0; 1000} 2000. . . cm2/vs" 
1450 input" be(min) cm2/vs " juuiud=<um-uu) /ny 
1460 printcb*jprintiprint" diagram drowing " 

1470 input" enter’ new data "jnd*:if nd$<>"n" and nd*<:>"d" then 1460 

1480 if nd* B "n" then 1520 
1490 if dH*="b" then 1510 

1500 gosub 1850iif dw* B "u" then 1520xrem gosub diagramm be=f(u(gs))a 
1510 gosub 2020irem diagramm sqr ( )=f(b> 

1520 gosub 1800: rem marker choice 

1530 if dw5="b" then 1640:rem onUaion of be a f (u (gs) ) -data 
1540 close6Jopen6,5jrem -data in be B f(u(gs>) 

1550 for gn =1 to gf 

1560 x=dx-int < (ga+(gn-l ) »gs) *dx/uw+,5> jy=int < <le4*b (gn, l)-uu) »dy/ud+. 5) 

1570 if y=0 then 1610 

1580 if >i<0 or y<0 then stop: gotol610 

1590 if x>n>stdx or y>ny*dy then stopjgotu 1610 

1600 print#6, "M"+str* (x+x 1 ) +" , "+str* <y+y2) :print#6,mn*- 

1610 next gn 

1620 print#6, "H";close6: if dw$="u"' then!790 

1630 printsprintjprint" sqr ( )«f(b) data "sgosub lSOOirem marker -choioe 

1640 cJo5e6jopen6 f 5:printch*:rem sqr < )«f(b> 

1650 pr int: print" choice of u (g* ) -value "jprint " B=e=nsne=3*===B=SHe=ac=s" 

1660 print:print"u(gs>= "jga}"bis " jge} "v" j spc (20) f " »te P width "}gs;" y" 

1670 printiinput" u(gs) =ca v";g 

1680 gn=int ( (g-ga) /gs+1+,5) :g=(gn-l) *gs+ga:print Jprint" u(gs)» "jgi"y" 

1690 print Jprint," wert ok": gosub 'HO; if j*<>"j" then 1650 
1700 for bn=2 to bf 

1710 x=int < (be-(bn-2) *bs) *dx/250+.5) :y=int (s(gn,bn) tny*dy/bm+. 5) 

1720 if x<0 or y<0 then stop: goto!760 

1730 if >:> nx*dx or y>ny*dy then stopxgoto 1760 

1740 if >:>0 and y~0 then 1760 

1750 print#6, "M"+str*(x+xl)+", "+str*(y+yl) sprint#6, mn$ 

1760 next bn 

1770 print#6, "H"sclose6 

1780 printch*;print:print" neues u(gs) ":gosub 930: j f j$="j" then 1630 
1790 close6:return 
ready. 
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1800 rem marker ♦(## 

1810 printsprint"markes "sprint" punkt ( 1 > "s print" rhombus (?)" 

1820 print' 1 square (3)":print" triangle (4) "sprint'double cone(5>" 

1830 print" octagon (6) "s input" ";mnSsmn*»''ir'i-mn* 

1840 return 

1850 clo6e6:open6, 5:rem ### be*= f(u(gs)) diagramm ### 

1B60 ua*< s "M"+str*(xl) + ", "+str*(y2) s printed, uo« 

1870 ux*="Xl, "♦str* (dx)+", "+str* <nx ) suy*="X0, *’+str*(dy) +" , "+str» <ny > 

1880 pr int#6,ux#s print 46, uyts print #6,uo*s pri ntM6, uy»:print#6,ux» 

1890 for n=0 to nx 

1900 print #6, "M"+str* <x 1+ntdx-hbt < 1+abs (sgn < 1 — n ) >)>■*■" , "+str< <y2-sd/2) 

1910 print #6, "P"+strt< (1-n) *uw> 

1920 next n 

1930 pri nt#6, "M"+str* <xl+nx*dx-10*hb)+", "+str* <y2-sd) 

1940 print#6, "S"+str* <sp ) : pr in t*)6, "P U(GS) "sprint#6, "S"+str*(sg) sprinted, "P/V" 
1950 for n=0 to ny 

1960 print#6, "M"+str*(xl-6*hb)+", "+,str$ <y2+ntdy-spt3) 

1970 pr int#6, *'P"+str* <uu+n*ud) 

1980 next n 

1990 print#6,"M"+str»(xl-7fhb)+", "+str« (y2+ny*dy*-10*hb> 

2000 print#6, H Ql"sprint«6 f "S"+str« <sp+l ) :print#6, "Pb"sprint#6, "S"+str«(sg) 

2010 print#6, "P/cm2/Vs“: print#6, "Q0" s pr i nt#6,uo«s cl ose6: return 
2020 cl ose6s open6, 5s rem ### sqr < ) K fib) diagramm ### 

2030 bo$'s"M"+5trt <x 1 )+", "+str*(yl ) sprint #6, bo* 

2040 bx*="Xl, "+str* (dx ) +" , " fstr* (nx) sby*“"X0, "+str* <dy ) +", "+str$ (ny) 

2050 print#6,bx*sprinttt6,by$:print#6,bo*sprint#6,by*sprint#6, bxS 
2060 for n=0 to nx 

2070 print#6. "M"+str$(int (xl+n*dx-.6»hb*ien(str«(n»250) ) ) )+", "+str« (yl-sd/2) 
20e0 print#6, "P"+5tr*(n*250) 

2090 next n 

2100 print #6, "M"+str* (:< 1 +n:s »dx-10*hb) +" , "+str* (yl-sd) 

2110 print#6,"S’ ! +strf <sp) sprint#6, "PB":print#6, "S"+str* <sg> :print#6. "P/mTesl a" 
2120 for n=0 to ny 

2130 print#6. "M"+str* (x l-6*hb) +" , "+5tr* (yl+n*dy-sp»3) 

2140 printM6, ”P"+str* (n<bm/ny ) 

2150 next n 

2160 rem root sign 

2170 x2=xl-7*hb:x3=int <l.B*hb+.5> :x4=int (2*hb+,5) 

2180 y3=yl+nyldy“12*hb:y4=(sg+l) *4 

2190 print#6, "M"+str* <x2-x3) / r", "+str« (y3) sprint #6, "10, ”+str*(y4) 

2200 print#6, " I "+str* (x3) +" ,0, "+5tr* <-2»x3> "+5tr$ <y4) 

2210 print#6, "10, "+ 5 tr«(9*hb)+", "+str9 <y4)+", 0" 

2220 print#6, "M"+str*<x2-x4> "+str* <y3+4»y4> sprint#6, "Ql" sprint #6, "PS(0) “ 
2230 print#6, "R"+str« (int (hb/2+.5) >+", 0"sprint#6, "10, "+str« (-4*hb) 

2240 print#6, "M"+str* (x2) +", "+str* <y3+4*y4) sprintttO, "PS(B> " 

2250 print#6, "M"+str$ (int <x2-x4/2+.5) ) +", "+str*(y3+4*y4+5*hb) :print#6, "P-l" 
2260 print«6, "Q0":print#6,bo*sclose6sreturn 
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100 
no 
120 
130 
140 
150 
* 160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
•320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
. 670 
680 
690 


rem 4## 

rem ### s = f uu (gs) ,u Tds) a tiff# 
rem #*# 

open 16, I6:chl=chr*(147) :printch*s print: print# 16, "bB";f2=0 


print" con«tnnU 11 1 pr i n t : pr i n t 
prints input" date ";da* 

prints input" to , t "spn« 

printsrl®!. ljpriin. ' load roniaUnco ,/ohm "jrl 

printsprint" u (as) "sprint" ■ — " 

input" starling value /V "jgas input" end value /v "jge 

input" 1 step width » /v " 5 gs: gf«int < (ge'-ga) /gs) + 1 
if gf<l then gosub 1230: goto 180 

i f (ga-. 7) * (ge-. 7) <0 or <ge+5. 5) * (ga+5. 5X0 then gosub 1290sgotol80 

printsprint" u(ds) " sprint " " 

input" atorting value /V "jdajinput" end value /v "jde 


input" .tap width /v ";dt:df=int ( (de-da) /dt) +2 

if df<l then gosub 1230s goto230 

if df>9 then gosub 1260s goto230 

if gf *df >200 then gosub 1260;gotol80 

if abs(da)>5 or abs(de)>5 then gosub 1290s goto 230 

printsprint" alternating voltage "sprint" " 

input" frequency /|ihZ "}f 

input" ampl i tuae/mv "ja 

printsprint" neauuring range "sprint" — — " 

input" sensi ti vi ty/uv ";se 

input" voltmeter dc / v " ; br : br=20000/br 

input" voltmeter ac /mv ars ar*20000/ar 

dim s(gf ,df ) sdim sd* (gf + 1 ) : ge=ge+. 001 *as: de=de+, OOUdt 

printchlS printspf intSprint" ■enaurenent running " 

ma-ti sgn=0sdn=0s rem neaaurencnt etnrt 
for d=da to de step at 

print:print:print"u (ds) u(gs) s'-sprint 

dn =dn+lsgn=0 
for g-ga to ge step gs 
gn =gn+l 

gosub 800s gosub 1050 
rem r*s acfc,dc-dmm 
ne>:t g 
newt d 

pr int#16, "bS" 

ms=(ti-ma) /60: mm=i nt, <ms/60> s m*=str* (mm) + " min"+str$ (int (ms) -60*mm> +" sec" 
printsprintsprint" end of neaaureaent/ expression” " 

gn=- 1 : s 1 *= "CCCCCCCCSC " : s2*=“CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC " 
f2=f2+l:if f 2 =2 then print#6: pr i nt#6s pri nt#6, , "b= 1000 mt":goto 630 
open6, 6 

pri nt#6, , " transoonductonce neoaurement "}da$ 

print#6, , " =====================": pri nt#6: print#6sprint#6 

print#6, , " teat "jpn$:print*6 

print#6 f lead real, >tance _ ";ri;" ohm" 

print#6, ," alternating voltage "ja;" mv / "sf;" khs " 

print#6, , "sensitivi ty "{sej" uv" 

pri ntD6, , " aeasuring tine . ";m< 

print#6:print#6, , " s=f (u<gs) ,u(ds> ) / ms N :print#6isd*<0)=" u(gs)2 " 
f l=3s fm=3:n<=-l:print#6, spc (24) j "u (ds) " 
for g=ga-gs to ge step gs 

gn=gn*i:if gn>0 then sd*(gn)=rightt <" "+str$(g)+" A ”,10) 

dn=0 

for d=da to de step dt 
dn =dn+l:if f 2=2 and gn=0 then 710 
if gr,=0 then sdt (gn) =sd$ (gn) +right* ( " 


+strf <d) + 


, 8 ) i goto ? 10 
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7C0 sd*(gn) a sd* (gn) +ri ght* < " "+str* (s (qn, dn) > + ” ‘‘,0) 

710 next 

7120 f l 3 f l+l : f m=f m+1 ( i f f 1 <5 then 760 

730 if fm=5 then s3t*l ef t* <sl*+s2*+s2*, len (Bd* (0) ) ) 

740 if f m< >5 then s3* a ri ght$ < ’’ A ’’,10) 

750 f 1=1 :print#6, so* 

760 printH6, sd* (gn) 

770 next 

780 if f2=l then print, nagnetic fieidjgosub 930: if j*=" j" then 380 
790 close6:printch*:priht:print:print:print, " ende ’'send 
800 rera )MMt 

810 rem MMM r!<s voltage source #)M) 

820 rem ### 

830 openl8,9 

840 ua*“"ual "+str* (g) i ub*="ubl “+str* (d) : print #18, ua* :print#18, ubt 
850 wz~60tgosub 1020 
860 if gntdnOl then 910 

870 i nputDIB ,r*trn=l en <ri ) i i f rn re 24 then rn B rn-l 
880 ua=val (right*(rt,rn-3) ) :ub e val (right* (r*,rn-15> ) 

890 if uaOg or ubOd then print ’’voitoge source on ";gosub 970: goto840 
900 print#16, "a8" 

910 clo5el8:printd,qj s return 

920 rem input loop 1 

930 get j*:get jSsget j*:get j$:get j*: j J="" 

940 print: print" ’ j/n' ?" 

950 get j* : if j$=’”’ then 950 
960 return 

970 rem input loop 2 

980 get j*:get j*:get j*:get j*:get j$: j * = H 

990 print:print" 'j' ?" 

1000 get j»: if then 1000 

1010 return 

1020 wt=ti irem wa!tlns 1o °p 
1030 if wt +wz > ti then 1030 
1040 return 

1050 rem MMM measuring value of ac & dc keithley dmm's MMM 
1060 open7,24:qpen20, 10: j=0; W2=300;gosub. 1020 
1070 input#20,w*; if stOOthen 1070 
1080 input#7, zi: if stOOthen 1080 

1090 j=j+l : i f w$="over"then print:print"overflow":6top 
1100 if z*=’’over"then print:print"overflow":stop 
1110 r4=r3:r3=r2:r2=rl :r 1-val <w*)/br 

1120 i 4=z3: z 3*22: z2=z 1 : z l=val ( 2 $) /ar : wz= 150; qosub 1020:if j<5 then 1070 
1130 if j !'• 1 0 then print:print" ideal value not convergent "lj“5:stop 

1140 r= (r l+r2+r3+r4) /4: rq=(r l-r ) 62+ (r2-r ) 62+ (r3-r) 62+<r4-r ) o2 
1150 z-(zl+z2+z3+:4)/4: zq= < 7 1-z) o2+ (z2-z ) 62+ <z3-z) o2+ <z4-z ) 52 
1160 if zq=0 then 1180 

1170 if z/sqr <zq/3><89 then 1070: rem 89: 95%wfc. for f< 17. 

1180 if rq=0 then 1210 
1190 ra=sqr <rq/3) srx=r/ra 

1200 if rx<89 then 1070; rem B9 : 957,w(t. for f< 17. (satz von laplace) 
1210 s (gn, dn) ~int 1 « OOtsofr/rl /z/10+. 5) / 100: print, 5 (gn, dn) 

1220 close?: cl osel'O: return 

1230 rem (MM) Sr error input 1 MMM 

1240 print:print" consult instructions "ipri nt:wz=60:gosub 010 
1250 return 

1260 rem MMM sr error input ,2 MMM , 

1270 print:print" anount „eaauring value too high "jprint: wz=60:gosub 810 
1280 return 

1290 rem ### sr °*-ror input 3 
lo00 print:print measuring range exceeded 

1 3 1 0 return 


f! :print:wz a 60:gosub 810 
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